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ABSTRACT 


Feedback techniques are employed to examine the influence of various 
physical design parameters on the aeroelastic stability of tilt proprotor 
aircraft. In addition, the influence of single loop feedbacks to improve 
the stability of the system are considered. Reduced order dynamic models 
are employed where appropriate to promote physical insight. The influence 
of fuselage freedom on the aeroelastic stability is examined as well as 
the influence of the airframe flexibility on the low frequency modes of 
motion relevant to the stability and control characteristics of the vehicle. 
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NOMENCLATURE 


•Conventional helicopter notation is followed in this report. Quantities 

are made dimensionless using air density, p, rotor rotational velocity, fl, and 

} 

rotor radius, R. The nondimensionalization is implied by the notation , e.g., 

V N 

the inflow ratio is written as V. Inertial properties normalized by 

are denoted by ( )*. Quantities are defined when introduced in the text or in 

Appendix I. The notation is identical to that of References 1 and 2. 



INTRODUCTION 


One of the critical problem areas related to the design and development 

of low disc loading tilting proprotor aircraft has been associated with 

eliminating the occurrance of aeroelastic instabilities at high cruise 

flight speeds. Owing to the importance of this problem area many studies 

1-12 

have been conducted related to this complex problem 

This report examines a number of aspects of this problem, the aero- 
elastic stability characteristics of tilting proprotor aircraft in cruise 
flight including the effects of fuselage motion. It has generally been 
found that a large number of degrees of freedom, as well as- precise modelling 
of the elastic properties of the rotor system are required in order to predict 

3 

all of the relevant aeroelastic phenomena . It has also been shown that these 

9 12 

models agree well with experiment ’ . The high order of these mddeis makes 

it difficult to obtain physical insight into the significant parameters in- 
fluencing the dynamic characteristics. 

One of the objectives of this study was to employ automatic control 
techniques to examine the feasibility of using relatively simple feedbacks 
to Improve the dynamic stability characteristics of the wing -proprotor system. 
It was also expected that through the use of these techniques it would be 
possible to obtain insight into the manner in which various physical param- 
eters influence the dynamics of this complex aeroelastic system. 

Many of the studies referred to above have been concerned with the 
dynamics of the wing prop-rotor system with the wing root assumed rigidly 
fixed, and therefore also of interest in this study was the influence of the 
fuselage degrees of freedom. The equations of motion of the vehicle were 
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extended to include the effects of the rigid fuselage degrees of freedom 

13 

in a previous study . The equations of motion for the complete vehicle 
system employed in this report are presented in References 2 and 13. 

A further objective of this study was to examine the effect of the 
flexibility of the wing proprotor on the dynamic stability character- 
istics of the vehicle as related to stability and control, and to examine 
the influence of body motion feedbacks as might be used for altering the 
stability and control characteristics of the aircraft on the aeroelastic 

stability of the vehicle. Experience on helicopters has indicated that 

14 

body motion feedbacks can influence rotor system stability 

The analytical model employed in the analysis presented is given 
by the equations of motion for the wing proprotor system presented in 
Reference 2. Appendix I presents the equations of motion taken from 
Reference 1. These equations of motion represent proprotor blade 
bending by an uncoupled model and are of simplex form than the model 
of Reference 2. The reduced degree of freedom models discussed in the 
first part of this report are based on the model of Appendix I. The 
refined equations of motion of Reference 2 taken with those of Reference 
13 are employed for the analysis of the fuselage free case. The modi- 
fications to the model of Reference 2, to account for body degrees of 
freedom are presented in Reference 13. 

Since one of the objectives of this study was to obtain physical 
insight into the dynamics of this vehicle, simplified dynamic models, 
based on Reference 1, are discussed in some detail prior to examining 
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the complete system. The reference physical parameters employed in this 
study are those of the Bell XV-15 presented in Reference 15. In some 
cases these parameters are varied from their reference values to obtain 

insight into the effects of changes in design parameters since other 

. 1 12 
proprotor aircraft designs have been considered ’ 

The topics covered in this report are as follows: First, the 

isolated rotor dynamics are discussed in cruise flight in order to 
obtain insight into the flapping and lagging dynamics of the proprotor. 
The discussion is largely concerned with the cyclic modes of motion. 
Cyclic proprotor flapping for the configuration of interest is a rigid 
blade motion owing to the gimbal mounting of the rotor, and cyclic 
lagging involves the first in-plane bending mode. Owing to the polar 
symmetry inherent in this problem, the method of complex coordinates is 
employed (Ref. 16) to reduce the order of the physical system from 
eighth to fourth order. Then the model is increased in order by allowing 
shaft flexibility. In order to obtain insight into this increased order 
system, first the dynamics of the proprotor on a flexible mounting is 
examined assuming that the support characteristics are isotropic. This 
allows retention of the method of complex coordinates and permits again 
a reduction in order of the system by a factor of two. The real physical 
system is then examined by relaxing the requirement that the support is 
isotropic. Various feedbacks are examined for the complete system 
indicating that certain decouplings exist with respect to the collective 
and cyclic motions of the rotor and their coupling to the wing motion. 

Then the fuselage degrees of freedom are relaxed and the dynamics of 
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the complete system with symmetric fuselage motion is considered in some 
detail. Consideration is also given to the use of reduced order dynamic 
models for the prediction of the motion of the vehicle relevant to its 
stability and control characteristics. 

A general arrangement drawing of the reference aircraft, the Bell 
XV- 15 is shown in Figure 1.5. 
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ISOLATED PROPROTOR DYNAMICS 


In this section the dynamic characteristics of an isolated proprotor 
are examined "with emphasis on high inflow conditions . Thus it is assumed 
that the support system of the proprotor is rigid and the dynamics of prop- 
rotor flapping and lagging motion are examined. The equations of motion 
development is presented in Appendix I. 

The proprotor physical characteristics employed for this study are 

3 15 

essentially those of the Bell XV-15 and are listed in Table I ’ . This 

2 

aircraft has a gimballed rotor with a soft flapping spring (v = 1.0355;) 

p 

so that its flapping motion is essentially that of an articulated rotor. 

The cyclic lag frequency, given in the Table, is the first flexible mode 
of the blade which in cruise flight is primarily lag bending (Ref. 3) . 

The frequency of this mode, as well as the relative proportion of flap 
and lag bending associated with the mode, is influenced by the trim air- 
speed as a result of the blade pitch variation required to trim the air- 
craft. As can be seen from Reference 3 it is essentially a lag mode at 
flight speeds above 300 knots. Therefore in the following it is referred 
to as a lag mode and the small component of flap bending is neglected. 

The complete equations of motion for the physical system assuming 
that the blade modes are uncoupled, taken from Reference 1, are given 
in Appendix I . The four degrees of freedom of interest in this section 
are the two components of gimbal motion or flapping, 3^^ and 3^^ and 
the two components of lag bending amplitude, and The equations 
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General: 


N = 3 


B == 12.5 ft 
f2 = 48 rad/sec 
=105 slug^ft^ 
Y =3-83 
a = .089 
a = 5.7 
h = .342 


Tor Cyclic Motion ; 

For Collective Motion 

- 

Constant RFM Case 

With RPM DOF 

vi = 1.0355 (*1) 

p 

V- = 1.85 

( - ) 

V = 1.33 (*2). 

’t? ^ ■ 

\ = ‘’ 

1$ = .670 

Jg = .779 
X-; = .670 

( -e- ) 

J* = 1.0 

= 1.035 


’•0 


(*1) : @ fl = 48 rad/sec 

C*2): @ F = .844 & P. = 48 rad/sec 


Table I. Description of the Rotor Used for. Numerical Examples. 
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of motion can be written as 


[A^] 


= 


where are given in Appendix I. 

Since the cruise flight condition is of interest the matrices 

and possess certain symmetry properties which can be employed to 

reduce the order of this system from eighth order to fourth order. A 

13 

set of complex coordinates are defined , 

5 ■ »lc " i * i 


S - Bjc i 

^IS 

c = ?. 

0=9- + i 

Ic 

®ls 

g = ct, 


IC XS g g 

Now by multiplying the second of equations [1) by i and adding to the 
first and multiplying the fourth equation by i and adding to the third 
we obtain two second-order equations with complex coefficients in place 
of four second-order equations with real coefficients. The equations of 


motion in this form are 


s - (Ma + 2i) s 

p 


+ (Vg - 1) + i 


- Qp (S - i) 


(s - i) 


-f 


s + (Q^ - 2i) s 


- (V- - 1) - i 


Mfl • iVM, 


% Kg 
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where 







The coordinates B and ^ are essentially complex multiblade coordinates 
and as such are expressed with respect to a fixed or non-rotating 
reference frame. The equations of motion can be effectively converted 
to a non-rotating frame and also incidentally, the complex coefficients 
removed, by defining a new operational variable 
i = s - i 

Substitution of this relationship into equation (.3) yields a simpler 
form of the equations of motion. 



where 

= - flj i * 

^ -2 ^ - 2 
CA J = s •+ Q- s + V 

^ uc ^ ^ 


These quadratic factors give the uncoupled flap and lag eigenvalues of the 

proprotor with respect to a rotating frame.- The eigenvalues are shown in 

Figure 1 as a function of advance ratio. Approximate expressions for the 

aerodynamic coefficients Mi and Q* are given in Reference 1, 

p Q 
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COS (j) 
8 


• 2 , 

sin (j) 

8 cos ({) 


where the effective radius approximation has been employed to obtain 
simple expressions for the aerodynamic derivatives. The inflow 
angle (j> is related to the advance ratio by 

V = - tan d) 

■ 4 ^ 

Thus from Figure 1, it can be seen that the flap damping, which 
is proportional to Mi decreases with increasing advance ratio and the 

p 

lag damping, Q* increases. 

As the model is complicated by more degrees of freedom it is con- 
venient to retain the original form of the equations given by equation 
(3) . The eigenvalues for this system are readily found from the re- 
lationships between the two operational variables s and s and thvis 
the fixed frame dynamics are found by adding i to the eigenvalues of 
equation (4) giving the values shown in Figure 2. 

.1.) Uncoupled Flapping Motion with Feedback 

First we examine the uncoupled flap dynamics in further detail 
especially as regards the response characteristics as a function of 
advance ratio and the influence of cyclic feedback on the flapping 
response. 

The uncoupled flapping equation of motion is given by 


{s - (Ml + 2i) s + i Ml} M Mo 5 + iVM, g 

P P D U 


(5) 
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Figure 2. Uncoupled Flap and Lag Characteristic Roots (Cyclic Modes 
Fixed Frame. 
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The influence of the weak flapping spring as represented has 

been neglected since its influence is small. The approximate expressions 
for Mg and are given by Reference 1 as 

M = 1 

8 cos t}> 

M = sin t{) 

U 6 

Thus the steady state flapping response to cyclic and gust inputs is 
given by 



( 6 ) 


In terms of the approximate aerodynamic derivatives, equation (6) can be 
expressed as 


^ss^ 0 - tan^4> g 

cos <})_ 

Thusi the flapping amplitude per unit cyclic increases with inflow angle 

_2 

proportional to cos <|) and the flapping response to gust inputs grows 
2 

as tan ^ . At an advance ratio of 1 (355 knots) using the 
equivalent radius approximation, the flapping amplitude per unit control 
deflection is 2.76 degrees per degree, and the flapping response due to 
a gust is 1.76 degrees per degree. The amplitude continues to grow 
rapidly with advance ratio as shown in Figure 3. 

Flapping feedbacks can be used to reduce the sensitivity of the 
rotor. Two possible flapping feedbacks considered are a 6^ hinge or 
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or pitch-flap coupling with the gearing denoted by K^, and an Oemichen 
type feedback where the pitch of one blade is controlled by the flap 
angle of an adjacent blade and the gearing is expressed as K^. In overall 
rotor plane sense, the Oemichen hinge is equivalent to a feedback which 


is 90 out of phase compared to the effect of the 6^ hinge. In complex 
notation, these two feedbacks can be expressed compactly as. 


0 = C- K + iK ) B C8) 

p o 

A root locus can be sketched as shown in Figure 4 showing the 
influence of these two feedbacks on the flapping dynamics. Of particu- 
lar Interest are the flapping stability limits and the amplitude reduction 
possible through the use of these feedbacks. Inserting the feedback law 
into equation (5j 

{s^ - (Mg + 2i) s + i (Mg - Mg K^) + Mg K^} B =.iVM^i (9) 


The influence on the flapping amplitude is given by the steady state 
solution to this equation 


g 


i V 


Ma 


i (1 - 


M„ M 

-®K j 


M 


3 


M 


B 


( 10 ) 


Substituting the equivalent radius approximation for the aerodynamic 
derivatives, the ratio of flapping to gust input amplitude is 
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{(1 + 



tan^ (j> 

K I - K 

° -)' * c 


COS (ft 


cos (f) 


Cil) 


This indicates that the amplitude reduction due to the 6^ effect 

(K ) is independent of the sign of the pitch-flap coupling and in addition 
P 

that both feedbacks become more effective as the advance ratio or 
inflow angle increases. 

The gearing ratios or feedback gains are limited by stability 
considerations. It can be shown that the stability boundaries 
determined from equation (9) are given by 

K K 

l + x cos(j>C — ^3^ - c— ^ 0 

cos <|) cos ^ 


Thus for the (S^ feedback only the stability boundary is given by 

K = - — • cos (j) 

P Y 

and for the Oemichen feedback alone 

K = cos^ <}> 
o 

Figure 5 shows the amplitude reduction which can be achieved by these 
two feedbacks, as well as the stability boundaries, indicating the 
effectiveness of these two feedbacks. 

The effectiveness of pitch-flap coupling in reducing the flapping 
amplitude can be seen, 

2.) Flap-Lag Coupling 

We now return to the complete set of flap-lag equations to examine 
the influence of flap lag coupling arising from the terms and Q^. 
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The coupled characteristic equation for the flap- lag motion can be 
written as in terms of the operational variable s 

uc , uc 

Root locus techniques can be employed to illustrate the influence of the 
coupling terms. Figure 6 shows the effect of this coupling term on the 
modes of motion- in the non-rotating frame. It can be seen that the 
coupling between the modes destabilizes the flap mode and stabilizes 
the lag mode. 

The proprotor is also equipped with a 6^ hinge- or pitch-flap coupling. 
This hinge causes additional coupling between the flap and lag motion 
resulting in a change in the eigenvalues of the coupled flap-lag motion 
as shown in Figure 7 at the reference advance ratio V = 0.844. It can - 
be seen that a negative value of (up flap causes and increase in pitch) 
stabilizes the lightly damped mode and the value for the vehicle (K^ = 

- 0.344) at this flight condition produces the maximum damping for the 
lightly damped flap mode. A further source of coupling arises from 
torsional flexibility of the blades. This effect can be approximately 

3 

treated as a pitch- lag coupling . For a more precise approach the coupled 
equations of Reference 2 are required. Figure 8 shows the effect of 
pitch- lag coupling on the eigenvalues indicating that the negative value 
(corresponding to a pitch reduction with increasing lag angle) reduces 
the damping of the lightly damped flap mode. The trend indicated is that 
expected from experience with helicopter rotors, although because of 
the fact that the flap mode has less damping it is this mode that is 
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Figure 6. Influence of Flap-Lag Coupling on Flap-Lag Dynamics 
i^Cyclic Modes) 
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Figure 7. Influence of Pitch-Flap Coupling on Flap-Lag Dynamics 
(Cyclic Modes) 
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Figure 8. Influence of Pitch-Lag Coupling on Flap-Lag Dynamics 
(Cyclic Modes) 
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destabilized by this coupling in contrast to the helicopter at low inflow 

where the lag damping is small and therefore it is the lag mode which is 
17 

destablized 

This completes the discussion of the isolated rotor cyclic dynamics . 

The collective flap and lag dynamics are also of interest. The collective 
motions can be treated separately from the cyclic motions . It should also 
be noted that the collective flapping is an elastic mode owing to the fact 
that the rotor is gimballed and each blade is not free to flap individually. 
3.) Collective Flap-Lag Motion 

The equations of motion describing the collective flap-lag dynamics 
are quite similar to the cyclic equations of motion expressed in terms of 
the modified operational variable s. The major difference arises from the 
fact that the collective flap mode is an elastic mode with a natural fre- 
quency, V„ '= 1.85. Pitch-flap and pitch-lag coupling are present. The 
effect of a gust disturbance enters in a somewhat different fashion. In 
addition, an additioijal degree of freedom, the rotor RPM must be included. 

In principle, the complete equation of motion for the rotor RPM degree of 

freedom would involve consideration of the drive train and engine governor 
19 

dynamics . It has been shown that these influences are not particularly 
significant and that the major effects of the rotor RPM degree of freedom 
can be examined through the study of two limiting cases: constant RPM, i.e., 
the case of a perfect engine/govemor system; and autorotation, where the 
rotor speed perturbation is determined by the aerodynamic forces acting on 
the rotor. 

With the constant RPM assumption, the elastic lag mode is excited. 
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however with the RPM as a degree of freedom this mode is not important. 
Consider the following two equations of motion, one describing the 


elastic lag motion. 



a ] + 
z 



AQ 


(13) 



Equation (13) reduces to, 

C _ a + vj ? = 4- AQ 

o z C ol* ^ 

• o o 

If the steady state torque is zero, equation (14) is, 

5 _ S = _i_ AQ (16) 

o z I* ^ 

0 

Comparison of these two equations of motion indicates that = 0 and 
therefore if the RPM degree of freedom (a ) is present, then the elastic 
mode degree of freedom (i;^) is not excited. Consequently the two cases 
considered are: RPM constant with elastic flapping and elastic lag, and 
RPM free with elastic flapping and no elastic lag. 

The equations of motion in the RPM fixed case are 
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uc 

- V 




. 

Co 


Mg VM^ 


% " 


(A D = s^ - M^ s -H v^ 
o uc ^o 

CA, ) = s^ + Q^ s + 

'“' ' «o 


C ) - I! 




Figure 9 shows the effect of the aerodynamic coupling as well as pitch-flap 
and pitch-lag coupling. The basic flap- lag coupling arising from aero- 
dynamics as well as the pitch-flap coupling act to reduce the damping of 
the collective flap mode. Pitch-lag coupling has only a small effect on the 
collective flap mode, and reduces the frequency of the collective lag mode. 

The free RPM case can be treated by taking v = 0 and considering 

^o 

the second motion variable as rather than C 3 i.e., - - t > Note 

s o s o 

that since there is no spring in the RPM equation the order of the system 
is reduced by 1. The uncoupled RPM mode has only one eigenvalue, s = - Q*. 

Figure 10 shows the effect of coupling in the RPM free case. Both 
in the case of the aerodynamic coupling and the pitch-flap coupling, the 
general trend is to reduce the natural frequency of the collective flap 
mode and produce only a minor change in the time constant associated with 
the RPM degree of freedom. Pitch-lag coupling is not present in the case 


with the RPM degree of freedom. 
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Poles are from Fig. (a) 


(s-plane) 



Poles are from Fig. (b) 



(s-plane) 



(a) Influence of Flap-Lag Coupling 


O; = -.233 


(b) Influence of Pitch-Flap 

Coupling; 

O: K - -.344 
P 


(c) Influence of Pitch-Lag Coupling 
(with Kp =5 -.344); 

O : K « -.3 


Figure 9. Collective Flap and Lag Motion Dynamics 


V = .844, Y = 3.83, 1* = . 779, 

^o 

I* = .67, v„ = 1.85, V = 1.33 














O : 166 


Figure 10. Collective Flap and 



O : Kp = -.344 . 

V = .844, y = 3.83, I* = . 779 

P_ 

I* = 1.0, Vg = 1.85, 


RPM Motion Dynamics 
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This section has considered the isolated proprotor dynamics which 
enter into the vehicle dynamics . In axial flight, the cyclic modes can 
be separated from the collective modes. 

In contrast to the helicopter which operates at low inflow, it 
should be noted that the cyclic flap mode has a lower damping than the 
cyclic lag mode. As the wing/support dynamics are introduced in the 
next section, recalling the characteristics of these isolated modes 
will be helpful in interpreting the changes which occur in the flexible 
wing case. 
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PROPROTOR WITH ISOTROPIC SUPPORT 


In order to promote physical insight into the aeroelastic stability 
of the wing/proprotor system, in this section the dynamic stability 
characteristics of the proprotor on a two degree of freedom isotropic 
flexible support is considered before examining the specific physical 
system under consideration. That is, for the actual aircraft the 
aeroelastic flexibility of the wing is modelled by torsion, chordwise 
bending and spanwise bending modes. In this section we only allow the 
support two elastic degrees of freedom and assvime that the inertial 
and stiffness properties are the same in both directions. The approxi- 
mate model permits the use of the complex coordinate transformation 
employed in the previous section to be extended to the case with a flexible 
support and helps to promote physical insight into the dynamic character- 
istics of this complex system. 

It should be noted that with this support model, the collective rotor 
degrees of freedom will not be involved in the problem, as they are primarily 
coupled to fore and aft translation of the proprotor shaft as shown later. 

Thus, in general, we are considering a system with two cyclic flap 
degrees of freedom, two cyclic lag degrees of freedom, and two directions 
of support bending. Only aerodynamic forces produced by the proprotor are 
included. The use of complex coordinates reduces this physical system 
from a twelfth order system to a sixth order system. 

Introducing the equality of the inertial, danT)ing, and spring 
characteristics of the support into the equations of motion given in 
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Appendix I, such that 


I* = I* = I* 

p py px 


c* = c* = c* 

X y 


C19) 


K* = K* = K* 

X y 

and introducing a complex coordinate describing shaft motion as 

d = a - i a 
y X 


The equations of motion for the three complex degrees of freedom can be 
-written as 


‘''b>uc 


I 

C s - i) 


1 


- Qg (s - i) 

i h ii^ (s - 1) 

2 

^ Cv^ - 1) 




S* 2 


- i h ) s 


- i h C s - i) 


s^ + [Ml + i (2 + h M )] s - i VM 

P M M 


S* 

i h s^ + (Q^ + i h Q^) s - i VQ^ 




(A 3 

a uc 


where 


'I 


a 


I* 

6 


1 and 


Mq iVM^ 

% 

- i h Ho hVH„ 
0 M 



C20) 
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M 


C ) 


ym 


( D 


^3 



YQ 


C ) 


I* 


I* 


I* + 2 Mg 




The tincoupled dynamics in each degree o£ freedom are 

* 2i) s (Vg - 1) + i 

4 CQj - 21) s 4 (v^ - 1) - i Qj C21) 

CA^) = A + (h^ H - i h Hx) s - h VH 

a-^uc p ^ p 3 y 

where 

2 2 
A = s + 2<;uj s + w 
P P P P 


with 



It can be readily seen that the equations of motion reduce to equations (4) 
when a = 0. 

First we examine the dynamics of the system with flapping and 
support degrees of freedom, then consider the influence of the lag degrees 
of freedom. 


Thus, eliminating the lag degrees of freedom from equations (20) 
the equations of motion are. 


t^3^uc 


I 

■H 


I - + (Mx + i (2 + hM )) s - i VM,. 

P PM 


i h (s - i) I 

" it 


(A ) 

^ a-^uc 


^6 
i h H 


iVM. 


0 


P 


h VH 


M 


3' 


a 


( 22 ) 
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The characteristic equation for this system can be expressed in terms of 
the characteristic equation for an articulated rotor system and a rigid 
propeller system as 

a = " ('’b - « <Vrp »35 


where (A„ ) is the characteristic equation for the articulated system 

3 + a fa 


^^6 + a^fa “ ^^3^uc 


(A ) 
a-^uc 


'3 = 1 


C24] 


+ i h (s - i) [s - (Ma + i C2 + hM )) s + i VM ] 
p p p U 


and CA^~) is the characteristic equation for the rigid system. 


I* 

(A ) = CA ) + .pi- [s^ - CMa + i (2 + hM )) s + i VM 3 

'• a-'rp ^ a^uc i* *■ g '• p-'-' p-* 


(25) 


It is interesting to note from the form of these two expressions that 
is the term responsible for coupling in the articulated case, and the 
ratio of blade to pylon inertia is responsible in the rigid propeller 
case. In both instances, the polynomial associated with the effects of 
coupling is the same and in a rather general sense can be- thought of as 
containing the effects responsible for the whirl flutter instability as 
will be shown. 

First, the simpler problem of the dynamic behavior of a rigid 
propeller described by equation (25) is examined. This leads 

to the classical problem of whirl flutter . I* may be taken as 1, and 
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the characteristic equation for the rigid case can be written as, 

= Ap + i (s^ + {hV^ - YM^ - 2i} s - ChVyH^- iVYM^))= 0 C26) 


It is convenient to use root locus techniques to illustrate the effect of 
the aerodynamic terms on the dynamics of the rigid propeller. By considering 
and fixed and varying the inertia, I*, we can examine the related 
importance of the aerodynamic te3rnis in causing whirl flutter. The poles 
of the system are given by A^ and with as the gain, the root loci shown 
in Figure 11 for three values of support frequency can be drawn. The pylon 
dynamics are shown, consisting of an advancing mode and a regressing mode. 

As the inertia is decreased, the roots move towards the complex zeros given 
by the quadratic factor in equation (26) . One of the complex zeros lies in 
the right half plane, moving farther to the right as the advance ratio is 
increased as shown in the diagram. The reduction in inertia (or conversely 
the increase is the aerodynamic terms) causes the regressing mode to become 
unstable thus, giving rise to whirl flutter. The advancing mode is stabilized 
by the coupling. Raising the support frequency also produces a stabilizing 
tendency. 

The zero in the right half plane lies there as a result of the term 
iVyM and thus it is possible to identify this term as the source of the 
whirl instability; it is interesting to note that this terra increases 
rapidly with inflow angle, i.e., from Reference 1, 

VM tan 4 sin i 

indicating that the onset of this instability tends to occur rapidly with 
increasing airspeed. As in more conventional flutter, the solution is to 
increase the stiffness of the mounting system, to reduce the effect of 
aerodynamics . 
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. Figure 11 . Root Locus : Influence of I* on Pylon Dynamics ; .Isotropic Pylon 
with Rigid Propeller; cf . Eqn. (26) 
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Now we turn to fully articulated case to see the contrast with the 


rigid case. Equation C24) is expressed as 
+ a^fa " = 1 ^ ^ 


(27) 


where 


j »» - ih-THj) s - WyH^} 

+ ihyHt (s - i) [s^ - {Mv + i (2 + hM ) } s + iVM ] 
p p U M 


(28) 


A root locus showing the influence of the support inertia on the modes of 
motion is shown in Figure 12. Again it may be noted that the critical case 
is the one with the lowest support frequency. The trend with Lock number, 

Y, is also shown, emphasizing the important role that the aerodynamic terms 
play in the instability. 

Now insight into the aerodynamic terms responsible for the instability 

can be seen by dropping the term from equation (28) . If = 0 in the 

fully articulated case, the blade flapping dynamics are uncoupled and only 

the support dynamics are influenced by inertia variation. Figure 13 shows 

this root locus. Comparison of Figures 12 and 13 indicate the important 

role of Hr in causing the whirl instability. It is also interesting to 
p 

note that in the articulated case, the instability occurs in both the advancing 
and regressing modes as shown by Figure 12. The location of one zero near 
the higher flap mode results in the fact that this root locus is essentially 
symmetric about the real axis giving a root locus diagram that is more con- 
ventional in appearance than in the rigid propeller case. That is, the 
influence of coupling is similar in both the advancing and regressing 
modes. It should be noted that Hi is of the opposite sign from its value 

p 

for conventional helicopter flight (low inflow) since it is equal to 
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Figure 12. Root Locus: Influence of I* on Coupled Dynamics; Isotropic Pylon 
with Gimballed Rotor; c£. Eqn. (27) . 
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Figure 13. Root Locus: Influence of I* on Uncoupled Pylon Dynamics (Hg 0) 
Isotropic Pylon; cf. Eqn. (28). 
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~ sin (j) , T 

% T“ ^ 

and therefore at high inflow conditions, the first term dominates. The 
negative sign corresponds to an unstable damping term at high inflow as 
comparison of Figures 12 and 13 indicates. With low support stiffness, 
the effect of can be seen to directly move the roots into the right 
half plane. Further insight into the fully articulated case can be 
obtained by considering an approximate factorization of This poly- 

nomial can be written as 

H, 


“fa = ' F » ‘ v|-» 

P = 1 u 


(293 


V 


Ma 


ii (s - i) {s - (1 - 


VM. 


-)} 




The simplified forms of the aerodynamic derivatives (Reference 1) give 

■ 1 


so that 


-“e 


’'“p 


3V 


HaM. 


= hV , {s - J (1 ^ {s^ - C2i (1 - ^ ) s 


y 


3V 


+ i M- n - — 


H Ma 

y B 


(30) 


Using the effective radius approximation 


^3^y ^16 

Vb ‘ 
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giving 


"fa ' - 1 - It "b’ - C21 * n V* 

The approximate zeros given by this approximate factorization agree well 
with the exact factorization of based on equation (28) . 

The last factor in equation (31) gives the zero very near the higher 
flap mode, the second factor gives the zero near the origin, and the first 
factor is responsible for the zero in the right half plane which may be 
viewed as the source of the instability, arising from Hx in effect. 

p 

Now we examine the influence of v^, the flap frequency, on the stability 
of the system. Figure 14 shows the influence of the flap frequency on the 
stability for three levels of pylon stiffness. Again with the high pylon 
stiffness, the dynamic system is stable and the effect of flap frequency is 
small. At a pylon frequency ratio of 1.5 the effect of is essentially 
to raise the blade flapping frequency as would be expected from the uncoupled 
system. In the low pylon stiffness case, the effect is markedly different, 
with the variation in causing primarily a damping change in the pylon and 
lower flap modes. In fact, it can be seen that there is an optimum flap 
frequency in the sense that the damping of the pylon mode and the damping 
of the lower flap mode can be maximized by a suitable choice of as 
suggested in Reference 7. For this numerical example the system is stable 
for 1.07 < v„ < 1.18. This figure also shows that the effect of the weak 

p 

flapping spring on the Bell rotor is small and that the dynamic character- 
istics of the system essentially correspond to the freely flapping rotor. 
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(b) 0) = 1.0 
P 



X :V3 = 

O :Vg = 

® :v^ = 


Figure 14. Root Locus: Influence of on Coupled 

Dynamics for Isotropic Pylon. 
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We now consider the influence of adding the lag degrees of freedom. 
The characteristic equation given by expanding equations (20) can be 
expressed in the following form 


3 + 5 + a 


' ^3 + ? ^P 




(32) 


where A . is the characteristic equation for the coupled flap-lag motion 
P 

and Ap is defined as before. If the inertia characteristic I* is infinitely 

large then the dynamics are given by A„ and A , the former associated 

3 + C P 

with the isolated rotor cyclic motion and the latter with the pylon. The 
root locus shown in Figure 15 presents the influence of pylon inertia on 
the dynamics of the system with the lag motion degree of freedom for the 
fully articulated case. The regressing modes are not particularly influenced 
by the coupling and the instability which occurs in the pylon advancing mode 
is quite similar to the case without the lag degree of freedom (Fig. 12) 
indicating that while the lag is required to obtain a detailed description 
of the dynamics, it does not play an essential role in the whirl flutter • 
instability. There also appears to be no particular tendency towards an 
air resonance instability undoubtedly as a result of the high lag damping. 
Figure 16 shows the effect of blade coiq)ling (blade pitch-flap and pitch-lag 
effects) which produces no essential change in the root locus diagram. 

Figures 12 and 14 appear to depict the essential contributing features 
to the instability. With a pylon frequency of 1.5, there appears to be 
little likelihood of an instability for the parameter ranges examined. 
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Figure 15. Root Locus: Influence of I* on Coupled Dynamics for Isotropic Pylon; 
With Cyclic Lag Degree o£ Freedom 
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1.) pylon Motion Feedback 

We now turn to the question of pylon motion feedback and its influence 
on the dynamic stability of the pylon proprotor system. First consider the 
rigid propeller case. The transfer function relating- complex pylon motion 
to complex cyclic pitch is found for the rigid propeller case by taking the 
limit Vg, ->- <» in equation (20) . The transfer function is 


a 

0 


YMq + ihyH 

I*iA } 
a 


rp 


0 


(33) 


where the denominator has been discussed previously. The feedback law of 

interest here is 

0 = - K.ct 
a 


The feedback gain can in general be taken to be a complex number to 
account for the phasing between pylon deflection and cyclic pitch. 

For example, if has a purely imaginary value so that = K^i then 
the feedback law implies 


01 =--K,a 
Ic d X 


01 = -K.a 
Is d y 


hence cyclic pitch is applied in such a way that the maximum decrease in 
cyclic occurs at the azimuth position (arg K^) greater than the position at 
which the maximum angular displacement of the pylon occurs. The character- 
istic equation for the closed loop system from equation (33) is. 
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This equation may be expanded to yield 


A 

P 


•J-t 


s(s - 2i) + 




U 


Mg) (s 


V(hH^ - iM^) + K^CMq + ihHg) 
h^H - Mi 

U 3 



C35) 


As seen previously, the possibility o£ whirl flutter is essentially due to 
the location of one of the complex zeros associated with the quadratic factor 
(see Fig. 11), i.e., which can be traced to the value of 

''n - " V 

l-X - Mg 

2 

and the sign of (h H - M^.) . Thus if the feedback gain K, is appropriately 

y p cl 

chosen, such that the quadratic expression has no zeros in the right half 
plane then the possibility of whirl flutter can be eliminated. 

Figure 17 shows the location of the zeros as a function of the feedback 
gain indicating that when = i, the first order factor has its root 
in the second quadrant and then the zeros of the quadratic factor are in 
the left half plane as shown for the parameter values appropriate to the 
flight condition. With this choice for the feedback gain, the system is 
stable for any value of the inertia parameter I* as shown in Figure 18 and 
thus the possibility of whirl flutter is eliminated. Figure 19 then shows 
the influence of the gain on the whirl flutter. It can be seen from the 
previous figure that the purely imaginary value of is of interest. It 
is interesting to note however that as the gain is increased, the regressing 
mode is stabilized, and the advancing mode is destablized and that a gain 
magnitude of less than 1 is desirable. 
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•Figure 17. Locus of Zeros: Influence of a ^ 0 Feedback; Isotropic pylon 
with Rigid Propeller; cf. Eqn. (35) 
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Root Locus : Influence of I* on Pylon Dynamics with a 0 Feedback 

(K, = i) ; Isotropic Pylon with Rigid Propeller; cf. Eqn. (35) 
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Figure 18. 






Figure 19, Root Locus: Influence of c; 0 Feedback on Eylon Dynamics; 
Isotropic Pylon with Rigid Propeller; cf. Eqn- (35) 
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Now consider the effect of the same feedback law on the gimballed rotor 
with no lag degree of freedom. The transfer function in this case obtained 


from equations (20) is 

(n“) . 

a ^ ©"^fa 


0 


^^0 + fa 


(36) 


The denominator is given by equation (24) and the numerator is obtained from 
equations (22) , 


<>fa 


ihyH 


j/{s2. [21, R. (1-^)1 a 


(37) 




. 1 


Using the effective radius approximation equation (37) becomes. 


The closed loop characteristic equation from (36) 

(Ao ). + Kj (nJ). = 0 

'• 0 + a^fa d ^ 0-^fa 


is 


(38) 


Figure 20 shows the root locus for this feedback for various values of gain. 
It is interesting to note that one of the zeros essentially acts to cancel 
the higher flap mode and the other zero lies quite close to the origin 
indicating that this feedback acts more like a damping term than a stiff- 
ness change. If is real, the feedback law is 


0 

0 


Ic 

Is 


K. 


a 

y 


K , a 
d X 


48 




-Figure 20. Root Locus: Influence of a ^ 0 Feedback on Coupled Motion 
Dynamics; Isotropic Pylon with QiaibMlled Rotor, 

i 

of. Eqn. (38) 
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This would appear at first glance to be a stiffness effect, however, examin- 
ation of the nature of the inplane force indicates that the dominant term 
in the inplane force is proportional to flapping velocity (at low frequency) 
and since cyclic pitch produces a flap angle this would imply effectively 
that a feedback proportional to shaft rate is produced by this displacement 
feedback. Shifting the phase of feedback by 90° produces more of a spring 
or stiffening effect as would be expected from the above discussion. It is 
also interesting to note that in general, a feedback phase angle which tends 
to stabilize one of the modes tends to destabilize the other mode, showing 
behavior quite analogous to the rigid propeller case (Fig. 19) . 

Similar behavior is shown in Figure 21 where the effect of various phase 
feedbacks is shown with different levels of pylon stiffness. In general, 
the phase of feedback which stabilizes or adds damping to one of the pylon 
modes, destabilizes or reduces the damping of the other mode. 

It is interesting to note the similarity of the results of Figure . 

21 to the case examined previously in which the coupling term in the 
pylon equation is dropped (H^ = 0) . In this case the closed loop character- 
istic equation is given by 


<■0 

a uc 


-K^ihHg = 0 


C39) 


The root loci based on equation C39) , i.e., for a dynamically uncoupled 
pylon is shown in Figure 22 for pylon frequencies = 1.0 and 1.5. It 
can be seen that the effects of feedback are very similar to the results 
shown in Figure 21. Of course in this case, the stabilization is not 
necessary since the increased stiffness has already stabilized the system. 
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Figure 21. Root Locus: Influcence of 
a 0 Feedback on Coupled 
Motion Dynamics; Isotropic 
Pylon with Gimballed Rotor 
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Now as a final topic in this section we examine the effect of the 
flap frequency in conjunction with feedback. The open loop transfer 
function is 


a 

e 


N 


C40j 


g + a 

The numerator obtained from equation (22) is expressed as. 






( 41 ) 


cc 

where CNg)£g^ is given by equation (37). Figure 23 shows the root locus 
for the zeros of the transfer function (equation (40)) as a function of 
Vg. The closed loop system equation is obtained from equation (40) as. 


■* a * ''e = “ 


( 42 ) 


It should be recalled that the system poles given by ^ ^ are also 
dependent upon v„. Root loci are shown in Figure 24 as a function of 

p 

gain for various values of flapping stiffness with a pylon frequency of 0.5 

The feedback influences the pylon and flapping dynamics. The phase of 

the feedback which stabilizes the motion shifts as the blade flap frequency 

is increased as would be expected due to the changing phase angle of the 

blade response to cyclic. It also may be noted that as the flap frequency 

is increased, the least stable mode becomes the flapping mode which is 

stabilized by feedback at the expense of the pylon mode damping. 

Figures 25 and 26 show the effect of feedback with increased pylon 

stiffness indicating that as shown before essentially the uncoupled 

result is obtained at co =1.0 and 1.5. 

P 
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Figure 23. Locus of Zeros: Influence of Vg Change; Isotropic Pylon; 
cf. Eqn. (41) 
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In Figure 27, the articulated rotor case is considered with a 
combination of displacement and rate feedback. 

With a suitable choice of rate and displacement gain the system can 
be effectively stabilized as shown in Figure 27(c) . 

The lag degree of freedom is added along with the appropriate pitch- 
flap and pitch-lag coupling and a root locus is shown in Figure 28. In 
this case a detailed discussion of the zero location is not possible. 
However, it can be seen that the primary effect of the additions to the 
model is to make the 90° phase feedback more effective in damping the 
unstable advancing mode and in addition to markedly reduce the effect 
of the feedback on the regressing mode. If the pylon stiffness is 
increased, the decoupling shown previously does not seem to be present 
as also shown in Figure 28. The feedbacks have little influence on the 
regressing modes and the phase of feedback which tends to stabilize the 
pylon advancing mode, destabilizes the higher lag mode. Thus, with the 
complete model it appears difficult to provide stabilization for the 
regressing mode and this result indicates the importance of including 
all relevant degrees of freedom. 

Figure 29 shows the influence of rate feedback and, displacement plus 
rate feedback for the complete system. It is interesting to note that 
rate feedback alone destabilizes the higher flap and lag modes and has 
little effect on the pylon mode. The combination of displacement and 
rate stabilizes the pylon advancing mode while destabilizing the higher 
flap mode. Thus, the addition of rate feedback does not appear to be 
useful in this case. 
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(a) Vg = 1.1 


(b) Vp = 1.2 


(c) Vp = 1.3 


Figure 24. Hoot Locus: Influence of a ->• 9 

Feedback on Coupled Motion Dynamics 
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Figure 28. Root Locus: Influence of ct 9 Feedback on Coupled Motion Dynamics. 
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DYNAMICS OF PROPROTOR AND CANTILEVER WING 


In this section we now examine the dynamics of the specific physical 
system of interest, that is, a cantilever wing and proprotor in cruise 
flight. The primary modification from the previous section is the support 
model which iS a non- isotropic support, i.e, a wing with three degrees of 
freedom:-, spanwise bending; chordwise bending; and torsion. First we consider 
the isolated wing dynamics. 

1.) Isolated Wing Dynamics 

The wing is modelled using only a single mode in each degree of freedom. 
The modal amplitudes are denoted by for spanwise bending, for chordwise 
bending and p for torsion. 

The torsion and spanwise bending are coupled as a result of the relative 
locations of the center of gravity and elastic axes. The chordwise degree 
of 'freedom is essentially uncoupled from the other two degrees of freedom. 

The coupled spanwise a ending -tors ion equations of motion are ejcpressed as 


uc 




.A 

^iP 

fA ) 

■-P uc 



0 


(43) 


The uncoupled characteristics are given by 

(A ) = (I* +1* + 2M? h^3 s^ + rc* 

p uc pw py D p 

(A ) = (I* + m* + 2M^ yl + s 

<ll uc P D tw ‘tw ■’ 


* 

2 


YC .) s + K* - hVvH 
PP P U 


(44) 


+ (C * + y“ YH. 


2 

"^tw 


.y 


- yC 




2N^-/yQ^) 


s + 


K* 
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and the coupling terms are 


A „ = (S* + 2M* hy hyH s 

pq^ w T) ^tw ^tw ' y 

Aq,p = (S* + 2M* h + (y^ hyH - yC 0 s - y VyH - yC 

^ w D ■' tv y q^^p-^ "^tw y qj^p 


The characteristic equation for this system is given by 


+ P 


= (A ) (A ) - A A =0 

uc p uc pq^ q^p 


(49) 


For the physical system under consideration the uncoupled quadratic factors 
have the following characteristics (at V = 0.844) 

(A^ )^c " (s + 0.035 + 0.413i)(s + 0.035 >- 0.413i) 

^Vuc " 1.10i)(s + 0.055 - l.lOi) 

It may be noted that the damping of these modes, which arises from aero- 
dynamics is very small. 

The coupling terms are 

A = 8.50 s (s + 0.074) 


pq 


qiP 


= 8.50 s (s + 0.599) (s - 0.599) 


The characteristic equation may be written in root locus form to examine 
the influence of coupling, 

1 _ 0.249 (sKs + 0.074) (s + 0.599) (s - 0.587) ^ ^ 

(s^ + .070s + 0.172) (s^ + .109s + 1.20) 

The root loci shown in Figure 30 illustrate the influence of coupling between 
the modes. Figure 30(a) shows the dynamics with all of the terms included 


65 



0° Locus 


Im 


0® Locus 


Ini 




Figure 30. Root Locus: Influence of ’q^ ^ p Motion Coupling on Coupled Dynamics 



and 30(b) shows the dynamics with the aerodynamic terms neglected, indicating 
that the coupling is produced by inertial effects arising from the center of 
gravity displacement from the elastic axis of the wing. The effect of coupling 
is to raise the p or torsion mode frequency and to produce only a small change 
in the spanjvise or mode from its uncoupled value. Examining the eigenvectors for 
these modes is also useful with respect to interpreting some of the results 
of the feedback analysis which follows. 

For the q^^ mode, 

E- - 1.26 
'*1 

and for the p mode 


5i 

p 


0.09 


-178°i 


Thus in the q^ mode, p and q^ act almost in phase while for the p mode,qj^ and 
p are very close to 180*^ out of phase. 

The chordwise mode is essentially uncoupled from the torsion, spanwise 
motion and the effect of aerodynamic terms is very small so that the chordwise 
motion degree of freedom motion is described by 


(A ) = [I * + M* + I* n| ^ + 2MJ (y.^, + h^ m?,)] s' 

q 2 uc qw p px tw Tj '"'tw ‘tw-'-’ 


(46) 


+ C* s + K* = 0 

^2 %■ 


Numerically 


(A ) = 110 (s + 0.012 + 0.677i) (s + 0.012 - 0.677i) 

q 2 uc 

Thus there are three frequencies associated with the wing motion, two associated 
with spanwise bending and torsion 

0 )^ = 1.10 , 0 )^ = 0.413 (uncoupled) 
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and one associated with chordwise bending 

w = 0.677 
^2 

The support system is, of course, no longer isotropic, however certain 
similarities can be noted with the isotropic support treated in the previous 
section as will be noted below. 

2.) Wing Pxoprotor Dynamics 

First we examine the influence of adding the wing degrees of freedom 
to the rotor cyclic motion degrees of freedom. Figure 31 shows the eigen- 
values calculated for various degrees of freedom as follows : 

4dof; This model includes only the rotor cyclic degrees of 
freedom ^ic ^Is^ ’ eigenvalues are 

those discussed previously. 

5dof; The wing torsion motion (p) is added to the rotor cyclic dof, 

6dof: The wing chordwise bending dof C<l2^ added to the 5dof 
model. 

7dof: The wing spanwise bending dof (q^) is added to the 6dof model. 

Figure 31 indicates that when the torsion degree of freedom is added the 
higher flap and lag modes are altered and there is little influence on the 
lower flap and lag modes. Moving to the 6dof case, i.e., adding the 
chordwise dof has essentially no effect on any of the Sdof eigenvalues. 
Similarly, for the 7dof case the primary influence on the 6dof eigenvalues 
is to raise the frequency of the p mode. This is of course, just the 
trend discussed in the case of the isolated wing. 

The complete system also involves the collective dynamics of the 
rotor and the rotor RPM as degrees of freedom. Figure 32 shows 
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(a) 4, 5, 6.& 7 DOF Models 


Figure 31. Comparison of System Eigenvalues for Several Rotor 
Support Models ((? V = . 844) 
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© : 2 dof or 

©: 7 dof (same as in Fig. (a)) 

X: 9 dof (7 dof + + ?o'(°r flg)) 





the two degree o£ freedom dynamics associated with these additional degrees 
of freedom for two cases: constant RPM including the collective flap and 
collective lag modes, and the collective flap plus rotor RPM degrees of 
freedom. As noted previously, with the rotor RPM degree of freedom, the 
collective lag mode is essentially not excited and therefore not included 
when RPM is included as a degree of freedom. 

Also shown in Figure 32 are the seven degree of freedom eigenvalues 
from Figure 31, The dynamics of the complete nine degree of freedom system, 
i.e., including the coupling between the two collective rotor degrees of 
freedom and the cyclic-wing motion system are also shown. It can be seen 
from this figure that there is only weak coupling between these collective 
dynamics and the 7dof. Ihe eigenvalues of the seven degree of freedom 
system correspond closely to the same modes in the nine degree of freedom 
case. 

Now we consider in more detail the influence of the support or wing 
degrees of freedom on the modes of motion and in particular the influence 
of anisotropy of the support as well as the nature of the coupling. First 
consider the effect of the addition of the torsion degree of freedom on 
the blade motion. The problem will be formulated as in the previous 
section such that the influence of varying the torsional inertia of the 
wing is employed to evaluate the nature of this coupling. The equations 
of motion for this system can be written as. 
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Where is the characteristic equation for the flap/lag motion, (A ) 

^ uc 

the uncoupled torsion motion. The coupling terms are given in Appendix 


is 

I. 


This various terms are given as 


hj - 1) . hYH^ 

T 

- ^8c s 

- hy s 

C = 


- hy s 

rp 

YQj s 

S* h s^ + hyH-s 


S* h + hyQ^s - yVQ^ 


( 48 ) 


and 


(A ) 

^ p uc 


(I* +1* + 2Mf h^) s^ 

pw py b ■’ 


+ CC* + h^ - C .) s + K* - hVyH 

p y pp p ' y 


(49) 


Pitch-flap and pitch-lag coupling are not included in the above matrices. 

They can be readily included by modifications to the matrices A^^^ and 

In order to evaluate the influence of the rotor support system dynamics, 

given by (A ) with the aerodynamic terms omitted (denoted by (A ) ) , the 
P uc ^ o ■ 

coupled characteristic equation given by equations (47) can be expressed 
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The first term in equation (50) represents the dynamics of the system when 
the support has infinitely large inertia associated with the torsion degree 
of freedom of the wing. In this case, the system dynamics are those of the 
isolated proprotor cyclic modes and the torsion of the wing without aero- 
dynamics. We now examine the influence of the torsion inertia on the 
dynamics, maintaining the natural frequency and the damping ratio 
constant, following a similar approach to that used for the discussion of 
the isotropic support. 

The root locus varying is shown in Figure 33. It can be seen that 
the influence of the torsional motion is primarily in the higher flap and 
lag modes along with the wing torsion mode. This locus then represents 
the effect of coupling, in the case where the support has one degree of 
freedom. 

We now examine the effect of additional wing flexibility arising 
from diordwise or motion. It is interesting to examine this case as 
a departure from the isotropic case by formulating the equations in the 
following way. Let 



The equations of motion can be written as 
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i I* = CO (corresponding roof .locus gain is 0) 

^o 

O : " = 2.66 ( " -.0703) 


Figure 33. Root Locus: Influence of I* Change on Coupled Motion Dynamics 

for 5 DOF Model (@ \T = .844; Z = .015 & to = 1.16) 

^o 
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It may be noted that if = (%>) , that is, if the support is isotropic, 

uc 

then the method of complex coordinates employed in the previous section could 

be applied. So it is convenient to examine the real or anisotropic case by 

splitting (Aqp_ into two parts, that is, 
uc 

CA ,) = (A ) + (AI* + AC*s + AK*) (52) 

where the quadratic factor in parentheses accounts for the anisotropy of the 
support. Thus the characteristic polynomial for the system, given by equation 
(51) can be expressed as 

^^3+ 5 + p = 

+ (AI* s^ + AC* s + AK*) (Ag + j. + p) 

where Ag ^ ^ p is the characteristic polynomial for the five degree of freedom 
case previously discussed (equation (47)) . The first term in equation (53) is 
the characteristic equation for the case of isotropic supports. Once the eigen- 
values are found for this system, then root locus techniques can be employed to 
examine the influence.;of the relevant inertia on the isotropic case and then the 
effect of anisotropy is introduced. 

The isotropic case is studied by introducing the complex coordinates 
as before with the addition of 
a' = p - i q^ 


11 
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rp 
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^ 2 ^ 
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hy s 
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The characteristic equation obtained from the first term in equation (53) takes 
the form of the previous section in the isotropic case, 

^6 ♦ C * C ■ * S 'V° ■ 'B '*? 

where ^ ^ is the isolated rotor cyclic flap-lag characteristic equation 
(12), and (4^)^ is the isotropic support characteristic without aerodynamics 
with the chordwise stiffness and inertial characteristic (qp the same as p . 
Figure 34 shows how the system dynamics vary with the support inertia 
again maintaining the natural frequency and damping ratio at their proper 
values. The inflow ratio V = 0.844 and the pitch- flap and pitch-lag coupling 
are included (K^ = - 0.344 and = - 0.3). For the proper value of the 
inertia, constant (I*^ = 2.66) the eigenvalues are shown. Note that since 
complex coordinates are employed here, the pole-zero configuration is not 
symmetric about the real axis. The trends shown are very similar to the 
isotropic results of the previous section as would be expected. Primarily, 
the advancing modes are influenced by the coupling and there is little 
influence on the regressing modes. The trend with inertia in the isotropic 
case is also very similar to the case with infinite stiffness in one 
direction as may be seen by comparison with Figure 33, since the rotor support 
regressing mode is hardly influenced by the coupling. 

Now the influence of the anisotropy is examined using the root locus 
approach described above. The factor under consideration related to the 
anisotropy is 

AI* s^ + AC* s + AK* = 33.8 (s + 0.021 +0.634i) (s + 0.021 - 0.623i) 


In the root locus shown in Figure 35, AI* is considered as the gain and 
AK* 

and are maintained constant. The root locus is based on equation (53) 
The root locus gain is given as AI* I*/ (I* I*^ - S* h) . Figure 35 then 
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Figure 34. Root Locus: Influ- 
ence of I* Change on Coupled 
P o 

Dynamics for Isotropic Rotor 
Support ((9 V = .844; =.Q15 

& u f= 1.161 ° 


* 5 ^ 1.01 


A 


Rotor Support Regressing Mode 


+-1.51 


(s-plane) 
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Figure 35. Root Locus; Influence of AI* Change on Coupled 
Motion Dynamics for 6 DOF Model (@ V = .844; 
AK*/AI* = C-623)^ & AC*/AI* = 2(.034) (.623) ) 
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shows the variation o£ the eigenvalues as a function of AI*. The isotropic 
case is represented of course by AI* = 0. The actual difference in support 
characteristics in the two directions is such that in fact the physical 
value of the inertia (AI* = 33.8) corresponds very closely to the limiting 
value of AI* = This case is the 5dof model and the chordwise parameters 
of the wing are such that the eigenvalues are close to this limit. 

The formulation presented is in general useful in studying the effects 
of support characteristics. 

This investigation appears to indicate the the wing chordwise degree 
of freedom for the parameters of the physical system of interest is not 
particularly influential. The chordwise mode does play a role however 
when the collective modes are examined as will be seen later in this section. 
3.) Wing Motion 'Cyclic Pitch Feedback 

We now turn to the examination of the influence of various single loop 
feedbacks on the dynamics of the wing proprotor system (nine degrees of 
freedom). An optimal control theory approach to this problem is presented 
in Reference 18. 

First consider the influence of q. 0, feedback at the inflow ratio 

^1 Ic 

V = 0.844. Positive (or negative) feedback represents a feedback law 
such that the increase in (positive for upward bending) results in an 
increase (or decrease) in 6^^^. That is, blade pitch takes a maximum value 
at = 0° (180°) in proportion to feedback gain. Figure 36(a) shows the 
root locus for the influence of this gain. Points are shown for a selected 
physically reasonable value of the gain. The gain can be interpreted as 
the ratio of the angular amplitude of cyclic pitch to angular deflection 
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(a) q-£ ©2c 


(b) 0^^ Feedback 


-. >HIh>-; Positive FB; : Negative FB, A&0: FB Gain = 2.51 (8 dB) 


Figure 36. Root Locus: Effect of 6^^^ or 0^^ Feedback (Proportional); 

9 DOF Model @ V = .844 


80 








o£ the wing tip measured from the wing root. It can be seen that the 
primary influence of cyclic feedback is to change the characteristics 
of the wing p and modes and the lower flap and lag modes. This feed- 
back essentially alters the damping of the torsion mode and the stiffness 
and damping of the spanwise bending mode. Positive feedback increases 
the torsion mode damping and destabilizes the spanwise mode. There is 
little influence on the wing chordwise mode, the higher flap and lag cyclic 
modes, the collective flap mode or the RPM mode. Considering a second 
feedback, wing bending to 0^^, i.e., a shift in phase of 90*^, essentially 
the root loci effects are shifted by 90°. In this case, the feedback 
alters the frequency of the p mode and the damping of the q^ mode. In 
each case, the larger change occurs in the lower frequency or q^ mode, and 
from the standpoint of stabilizing the q^^ mode, 6^^ feedback appear to be 
the best candidate. Now consider the effects of torsion motion feedback 
shown in Figure 37. The dominant effects of this feedback are to change 
the torsion mode characteristics as would be expected and there is some 
influence on the higher lag mode. A positive feedback of 6^^^ can be 
employed to provide torsion mode damping. It is interesting to note 
that the trends with torsion feedback are very similar to the influence 
of spanwise bending feedback with respect to their influence on the wing 
modes. Of course, the q^ feedback causes a larger effect on the q^ mode 
and similarly the p feedback causes a larger effect on the p mode. The 
other difference to be noted is that the torsion mode feedback produces 
favorable effects on both modes (although the effect on is small) 
whereas the q^ feedback produces opposite effects. This can be explained 
by noting the eigenvectors for the p, q^ coupled modes given earlier 
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X— Positive FB, • Negative FB, A& 0 : FB Gain = 1.0 (0 dB) . 


Figure 37. Root Locus: Effect of p -»■ 9^^ or 9^^ Feedback (Proportional); 
9 DOF Model @ V = .844 
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where in the mode, p and are in phase, whereas in the p mode they are 
out of phase. 

Again the p mode feedback has little influence on many of the modes 
similar to the q^ feedback. 

Torsion feedback is essentially similar to a single axis version of the 
focussed pylon mount. 

Now we examine chordwise motion feedback. Figure 38 shows the effect of 
this feedback to each cyclic. The effect is generally small due to the large 
stiffness noted above. For a reasonable physical value of gain some slight 
increment in the damping of the chordwise mode can be provided. 

These results indicate that owing to the small influence of q^ feedback 
a focussed mount system with equal gains on both axes would act essentially 
like a single axis feedback of torsion or spanwise bending motion. 

In considering more generally the influence of these feedbacks it is 
possible to note some marked similarities between the various feedbacks. 
Further it is possible to associate certain of the zeros with the wing 
dynamics and others with the blade^ dynamics. 

Figure 39 shows the zeros for the q^^ transfer functions indicating 
the source as arising from wing or rotor characteristics. The zeros from 
the seven and nine degree of freedom models are shown to indicate that there 
is only a small difference in these two cases. By the symbols H and W the 
origin of these zeros from the wing or H force transfer function are 


identified. Figures 40 and 41 show the zeros of the other transfer functions. 

"^1 ’"12 p 

Each of the transfer functions — , s — has nearly in common the 

®1C hs 'he 

following zeros (denoted by H , H ,Y , or Y ) : 

S C 3 O 
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(a) q- 0- Feedback 
^2 Ic 
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EPM Mode 
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Mode ^ 


Zero 

@ s = .74 


(h) q„ 6, Feedback 
2 Is 


X — E}-> : Positive FR^ : Negative FB, A&El: FB Gain — 2.51 (8 dB) 


Figure 38. Root Locus: Effect of q 2 -»■ 9^^ of Feedback (Proportional); 


9 DOF Model @ V = .844 
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(a) Zeros of (b) Zeros of q^/Sig 



Figure 39. Zeros of Transfer Functions, q^/^^c ^ ® ^ “ .844 


85 








(a) Zeros of p/0 


ic 


(b) Zeros of p/0 


Is 


W: Zeros Originating from Wing Dynamics 


H or H : 
c s 


Rotor 


(H Force to 0^^ or 6^^^) 


Figiire 40 » Zeros of Transfer Functions, p/^^c ^ ® ^ “ .844 
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(a) Zeros of q2'^®lc Zeros of 

W: Zeros Originating from Wing Dynamics 

Y or Y : " Rotor " (Y Force to 0, or 6, ) 

c s • ^ Ic Is 


Figure 41.' Zeros of Transfer Functions, ^2^®lc ^ ® ^ ~ 
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1.0 Three zeros on the real axis; one on the negative real axis 
near the origin Cs = - 0 . 1 ); another farther from the origin 
(s. = - 1,4); and one on the positive real axis (s = 0.8). 

2. ) One zero near the lower mode flap eigenvalue (s = - 0.1 + 0.14i). 

3. ) One zero near the point s = 2i. 

*^1 ^2 p 

Similarily for the other three transfer functions, - 5 — y - 5 — , there are 

0, On 0 

Is Ic Is 

1. ) Two zeros on the negative real axis; one near the origin 

Cs = - 0.05) and one farther to the left [s = - 0,5), 

2. ) One zero near the lower flap mode eigenvalue (s = - 0.03 + 0.2i). 

3. ) One zero near the point s = 2i. 

4. ) One zero far from the origin (s = - 0.6 + 5.3i). 

As will be shown below, these zeros in both cases are associated with the 
rotor dynamics, and are essentially the zeros of the rotor inplane force 
to cyclic transfer functions. The following section considers this 
question in more detail. 

4.) Transfer Function Zeros for Cyclic Feedback 

We now examine some characteristics of the zeros of the transfer 
functions of the coupled proprotor cyclic and wing motion to illustrate 
the source of the zeros. The location of these zeros governs the manner 
in which feedback influences the dynamics of the system. 


The coupled system equations of motion may be written as 

; C > C t C " 

11 ' rq 4 rq i rp 

I _ ^ 1 t 
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where 
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roA 
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Only the seven degree of freedom model has been employed, since it has 
been shown that the collective rotor motion and the RPM degrees of freedom 
do not influence the dynamics of this system. 

Writing the equations (56) in more compact form 



(57) 


and defining a new rotor motion variable 


x' = X + A* 
r r 11 21 


w 


(58) 


The equations of motion can be written as 



(59) 


'll 


= B 


11 


(60) 
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Now in this form,x^ may be regarded as the isolated rotor motion. Thus to 
find the input terms on the right hand side of the equation above, in 
particular as regards the zeros, they can be found from the isolated 
rotor transfer functions ; 

For Vg = 1, which is a good approximation to the proprotor under con- 
sideration and noting that 

>^tw ' 

C ^ = C 
q^^r pr h 


We can express the right hand side of equation (59) as 
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The left hand side of equation (59) can be sin^lified based on the fact 

noted earlier that, wing chordwise bending almost no effect on 

the eigenvalues of the rotor wing dynamics, therefore we may drop coupling 

terms such as C . Further q, motion only effects wing torsion so that 
rq., ^1 
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The hub moment contribution can be neglected due to the low flapping 
stiffness. 

As a result of these approximations equation (59) can be expressed 


as 
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Now the transfer function g— may be written as 

Ic 


qi q^ f % fx 
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From equation (64) the numerator of is 
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and the ntunerator 
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(67) 
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Therefore it can be seen that the 3 — transfer function zeros are given 

®lc 

by the transfer function and the zeros of the transfer function — . 


Ic 


Similarly the transfer function 


q f 
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q, f 

1^ ^ X 
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~ 

since ^ ~ 0 the zeros of g — consist of the same wing zeros as for the 
X Is 

•g — with the rotor zero contribution given by . 

Ic Is 


The approximate zeros given by this approach are 
From 


(A ) = 111 5 + 5.17s +49.5 

^2 


s =ZQ. 02 ± 0.67i 


From 


fA 3 - A = 0.481 + 0.265 s + 3.96 

p-'uc y q p 
tw -L 


s = - 0.28 ± 2.86i 


y 

The zeros of the transfer function 

®lc 


s = - 0.11, 0.77, - 1.37, - 0.11 ± 0.14i, - 0.01 ± 2.03i 


The zeros in the figure are labelled W or H as to whether they arise from 
the wing dynamics or the rotor dynamics on the figures. The physical 
origin of the zeros can be seen. The wing motion response to cyclic can 
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be considered to be a product o£ cascaded transfer functions, i.e, the 
product of the wing response to rotor force times the rotor force to 
cyclic where the zeros can be calculated based on a decoupling assumption. 


Now consider the torsion transfer functions 
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From equation (64) , we obtain the numerator of 
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Again the zeros are given by the wing contributions 
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and the zeros of the transfer function . Similarly 
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So the wing contributions are the same as above and in. addition we have 


f 

the zeros of the transfer function . 

®ls 

Thus the torsion to cyclic transfer functions have the sets of zeros 
from the rotor which are identical and differ only in the wing contributions. 
In a similar fashion we can examine the chordwise transfer functions 
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and the numerator of 
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so it can be seen that the zeros for chordwise feedback are: 


The zeros of the factor [ (A ) (A ) - A A ] which are in 

uc p"uc pq^ q^p 

fact the zeros of the characteristic polynomial A , i-e, the coupled 

qi + P f 

.X 

wing torsion bending motion. And the zeros of the transfer function - 5 — , 

f Ic 


which because of symmetry is equal to - . 

q„ Is 

Similarly the zeros for the transfer function have the same wing 

Is 

contributions (the uncoupled wing torsion bending zeros) and the zeros of 

f f 

the transfer function which is equivalent to The zeros obtained 

Is Ic 

from this approximation agree well with the results obtained from the 


complete model. 


Thus the source of the similarity in the zeros can be seen. To 


summarize the results. 


^1 ^2 p 

- 5 — ^ — all have essentially the same zero 

b, b, b, 

Ic Is Ic 

contributions from the inplane response characteristics and differ in the 

^1 ^2 p 

wing zeros. - 5 — , ^ — , and ^ — also have similar characteristics, indi- 

b- b, b. 

Is Ic Is 

eating that the zeros may be derived from a cascaded system in which the 
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coupling terms appearing in the numerators can be neglected, i.e., the 
zeros in the rotor force to cyclic transfer functions can be estimated 
ignoring wing motion. 

It should be noted that these conclusions depend upon a nvanber of 
assumptions which may not be true for other proprotor aircraft with 
markedly different characteristics. It appears that = 1 is one of 
the crucial assumptions . 

In this particular case, the difference in the action of each of these 
feedbacks is essentially due to the role of the wing dynamics. Chordwise 
feedback has little effect on the other wing modes owing to the fact that 
the approximate zeros of the chordwise transfer functions lie on the qj^,p 
modes . The torsion and spanwise feedbacks have one pair of zeros very 
near the mode. The other pair is basically related to the wing response 
to rotor force, and produces the nature of the feedback effect. 

To summarize, the essential difference in the effects of the various 
wing motion to cyclic feedbacks appears to be due to the location of the 
zeros due to the wing dynamics, and these can be estimated from the wing 
motion transfer functions in response to forces applied at the rotor hub. 
Both the and p feedbacks give a zero at the uncoupled chordwise mode. 

The second pair of wing zeros for q^ feedback is at rather high frequency 
while for p feedback the zeros are at relatively low frequency near the q^ 
mode, sufficiently close in fact such that p feedback has only a relatively 
weak influence on the q^ mode, q^ feedback results in zeros from the wing 
dynamics which are near the q^^ and p modes as would be expected. 

It is also interesting to note that the rotor force to cyclic transfer 
functions yield a zero quite close to the origin implying that the response 
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of the rotor force to cyclic is close to being a rate dependent effect. 

A simple model for the rotor in-plane force injilies that for low fre- 
quency motions the inplane force is proportional to flapping rate, and 
since cyclic produces a flap angle for low frequency inputs, the rotor 
force transfer function is more like a damping feedback than a stiffness 
change. This would be quite different if the rotor produced appreciable 
hub moments through a flapping spring. 

Thus it can be seen that cyclic pitch feedback could be effective in 
controlling the wing torsion and spanwise bending modes but has little 
influence on the chordwise mode. Therefore it seems desirable to examine 
the effect of chordwise motion to collective pitch feedback, 

5.) Chordwise Motion to Rotor Collective Feedback 

The root loci shown in Figure 42 illustrates the influence of pro- 
portional feedback of wing q^ motion to collective pitch for the nine degree 
of freedom model in the free rotor RPM case at advance ratio V = 0.844. The 
sign convention for positive Cor negative) feedback is such that an increase 
in <^2 (positive for aft bending) results in an increase (or decrease) in 
blade collective pitch. 

It can be seen from this complete model that essentially- only the 
chordwise mode, collective flap and RPM modes are involved and pole zero 
cancellations eliininate the effects on other modes. Thus it appears that 
this feedback can be adequately examined using only a three degree of 
freedom model. 

Figure 43 shows some results using this three degree of freedom model. 
It is interesting to note from Figure 42 the zero configuration involved 
in this feedback. There is a zero at the origin implying that this feedback 
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Figure 42. Root Locus: Effect of 9^ Feedback (Proportional); 9 DOF 

Model (in Free Rotor RPM Case) @ V — .844 
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is essentially a rate feedback. That is, that a collective pitch change gives- 
no steady state thrust. This is a consequence o£ the free RPM assumption. 

That is, in the steady state no thrust change is produced by a collective 
pitch change. The two real axis zeros indicate the presence of a non-minimum . 
phase effect. These trends can be readily explained by considering the three 
degree of freedom dynamics in more detail 

The equations of motion can be written as 
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Since this is a relatively simple system, we do not consider the cascade 
approach used in a previous section. 

The numerator of the -g— transfer functions, denoted by Ng is given by 
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But since the Q derivatives are proportional to the T derivatives^, i.e.. 
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equation (77) can be simplified to 
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Expanding 
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The damping term in the second factor is very small, i.e., if the effective 
radius approximation is employed this term vanishes, so that the numerator 
is given by 
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Inserting numerical values we find that the zeros are 
s = Oj + 6.49 

which compares well with the more exact result obtained by retaining all 
of the aerodynamic terras 

s = 0, - 6.06, 6.95 

The zero located at the origin gives a zero static gain indicating that 
a step input in collective produces no steady state chordwise deflection 
implying that in the steady state condition there is no thrust change. 

It can also be seen that the nonminimum phase effect, or zero in the 
right half plane arises from the blade inertia. A step increase in pitch 
input causes the blades to accelerate forward and the inertial reaction 
of the wing occurs in the opposite direction or to the rear. 

Consider the equations of motion for the three degrees of freedom, 
neglecting aerodynamic terms on the left hand side. 
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Figure 43 shows the influence of coupling for this model, including aero- 


dynamics, with and without the degree of freedom. The primary effect 

of coupling is to raise the frequency of the collective flap mode, as 

the neglect of the aerodynamic terms decouples the degree of freedom. 

It also can be seen that the S* term representing the blade inertial 

^o 

characteristics is responsible for the non-minimum phase term. 

The complete system analysis of Figure 42 shows an influence on the 
mode, a positive feedback, i.e., aft bending of the wing producing an 
increase in collective pitch results in a stabilizing effect. Figure 44 
considers the influence of rate feedback as well as a combination of rate 
and displacement indicating that a combination of rate plus displacement 
is most effective in increasing the damping of the wing chordwise mode 
although the damping of the collective blade mode is slightly reduced. A 
physically reasonable value of gain produces a rather small effect however, 
owing to the high chordwise stiffness of the wing. 

To summarize the results of this section, it generally appears that 
spanwise or torsion bending to cyclic feedback could be effective in 
controlling the dynamic characteristics of these modes. Owing to the high 
stiffness of the chordwise mode, it appears difficult to obtain much change 
in the chordwise mode owing to the high stiffness . Collective pitch feed- 
back appears more effective, for reasonable gain values, in controlling 
the chordwise mode. 

Other possible feedbacks have not been considered since they primarily 
influence the isolated modes, i.e., cyclic feedback to blade motion 
alters the blade motion itself and wing flap feedback will largely only 
influence the wing torsion-spanwlse bending motion. 
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Figure 44., Root Locus; Effect of 0^ Feedback; 9 DOF Model (in Free Rotor RPM Case) @ V = ,844 






PROPROTOR DYNAMICS WITH FUSELAGE FREEDOM 


In this section we examine the influence of freeing the fuselage on the 
overall dynamics of the proprotor aircraft with flexibility. Only the 
symmetric degrees of freedom of the fuselage are considered, i.e., pitch, 
vertical translation. and horizontal translation. Further, additional 
blade degrees of freedom have been included to more precisely represent 
the influence of the blade .torsion degree of freedom. The proprotor wing 
model is now that of Reference 2. Rather than using an effective pitch-lag 
coupling as in the previous sections, additional blade degrees of freedom 
including modal coupling between flap and lag bending as well as control 
system flexibility are incorporated. This increases the complexity of ' 
the basic model with free proprotor RPM. Fifteen degrees of freedom are 
now involved. Therefore with the fuselage free, we have an eighteen degree 
of freedom system. The inclusion of additional blade degrees of freeom, 
while changing quantitative aspects of the results does not give any 
essential qualitative differences from the nine degree of freedom model 
with pitch- lag coupling. 

Of particular interest here is the influence of freeing the fuselage 
on the dynamics and stability of the overall system especially as regards 
its influence on the whirl flutter speed. In addition, the influence of 
the flexibility on the dynamics of the aircraft, especially with respect 
to the lowest frequency modes which are related to the stability and control 
of the aircraft are examined and an approximate model discussed which con- 
sideres the primary influence of flexibility. 

Figure 45 shows the eigenvalues of the eighteen degree of freedom 
Cfuselage free) and fifteen degree of freedom Cfuselage fixed) motions 
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6 High Frequency Mode 
Eigenvalues 


6 High Frequency Mode 
Eigenvalues 
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Figure 45. Root Locus for Velocity Sweep; 
@ Free Rotor RPM 
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(b) @ Constant Rotor RPM (17 DOF Model) 

Figure 45 . Locus for Velocity Sweep; Aircraft Longitudinal Rigid Body Motion Modes 
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(a) 17 DOF Model; With Aircraft Longitu- (b) 14 DOF Model; Without Aircraft Longi- 

dinal Rigid Body Motion DOF tudinal Rigid Body Motion DOF 

Figure 45. Root Locus for Velocity Sweep; 

@ Constant Rotor RPM 
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as a function of trimmed flight velocity and also the case with rotor RPM 
assumed constant (seventeen and fourteen degrees of freedom) . It can be 
seen that there is little influence on many of the eigenvalues with the 
exception of the wing bending modes and of course the addition of the class- 
ical longitudinal modes of motion of the airframe. 


In particular it can be noted that upon freeing the fuselage both of 

the wing bending modes (q and q ) are raised in frequency and their 

^1 ^2 

damping ratios are increased. This increase in frequency is simply 
explained by considering the way in which the natural frequencies of a beam 
undergoing symmetric vibration are influenced by a mass located at the 
center of the beam. If we compare the lowest natural frequency of a 
cantilever beam which is given by 



which would correspond to the limiting case of an infinite mass at the 
center of a free-free beam to the natural frequency of a free-free beam 
with no mass at the center, its natural frequency is 





El 

4 

mr 


indicating the lowest mode natural frequency increases as the effect of 
fuselage freedom is included. This effect can be seen in the complete 
model by comparing the frequencies of the q^ and q^ modes in Figures 45 (a) 

Vi ^2 

and (b) . Figure 46 shows these results more explicitly giving the frequency 
of there two modes obtained from the fifteen and eighteen degree of freedom 
models. 
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Damped' Natural Frequency 
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Also shown is the influence of freeing the fuselage on the damping of the 
two modes. There is little effect shown on the damping of the chordwise mode 
however the damping of the spanwise mode is increased quite significantly 
apparently obtaining some additional damping from the freeing of the fuselage 

The other case of interest here is that with the rotor RPM constant, i.e., 
the case of perfect engine governor. The eigenvalues in this case with the 
fuselage fixed and fuselage free are also shown in Figure 45. The general 
trend of the eigenvalues with speed are similar to the RPM free case. There 
is a difference in the collective mode for reasons explained earlier. 

As far as the stability of the system, or in other words the whirl flutter 
speed, there is little change in the damping of the mode as shown in 

Figure 46. There is however a marked increase in damping of the q^, mode. 

1 

This increase in damping arises apparently from the fact that wing spanwise 
bending causes a rotation of the rotor shaft with respect to space and thus 
if the proprotor RPM is assumed constant aerodynamic damping appears in the 
spanwise mode. If the RPM is free, then in effect no aerodynamic forces 
appear as a result of this rotation. 

The constant RPM assumption also has a marked influence on the proprotor 
contribution to the phugoid damping as may be seen from Figure 47. At the 
higher flight speeds the phugoid becomes critically damped due to the large 
thrust variation with airspeed which does not appear for the free rotor RPM 
case. 

Figure 47 compares the time to half amplitude in these two cases for the 
phugoid and the short period where it can be seen that the only significant 
difference between the eighteen and seventeen degree of freedom models 
occurs in the phugoid mode. The period change in the phugoid between the 
two cases is of course due to the difference in damping. 
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Figure 47. Phu'goid Mode Characteristics 






"QS" Indicates that the Quasi-Static Approximation is Applied to All Motion DOF except the 
Aircraft Longitudinal Rigid Body Motion DOF 


Figure 47. 


Short Period Mode Characteristics 





We also examine at this point the influence of the aeroelastic degrees 
of freedom on the two low frequency modes associated with the stability and 
control characteristics of the aircraft. Also shown in Figure 47 are 


the quasistatic approximations to the low frequency motions, i.e., the 
dynamic terms associated with the aeroelastic effects are eliminated. 

This approximation must be made with care since we must retain the proper 
rate dependent aerodynamic contributions from the various modes. It can 
be seen that the only significant difference between the quasistatic 
approximation and the complete system dynamics is associated with the 
short period damping. Use of the quasistatic approximation results in an 
overestimation of the short period damping or time to half amplitude. 

It would be highly desirable in order to obtain insight into the aircraft 
dynamics to obtain a low order approximation to the short period motion 
including the effect of aeroelasticity. The following section considers 
this question in detail. 

1.) Short Period Motion 

First we consider the quasi-static approximation to the short period 
motion in order to justify an approximate way in which to account for the 
dynamics of the aircraft associated with its flexibility and then proceed 
to develop this approximation. 

The short period equations of motion in an inertial frame of reference 


may be written as 


Cl - Z^) J - 2^ 


- M. s - M 
w w 


- s - Zq 
- %s - Mq 



1 9Z 

( ) ~ M 3( ) 



1 3M 

I 9C ) 




where 



These equations are written in dimensional real time (s = sfi) . It is of 
course assumed that the center of gravity of the aircraft is fixed at 
one point and thus M is the total aircraft mass and the moment of inertia 
includes the rotor contributions as a point mass. Thus the rotor inplane 
force does not contain any inertia terms. The rotor blade inertia terms 
are included in the moment. . This is important to keep in mind 
when we examine the dynamic form of the inplane force and moment in the 
following. 

Note again that the conventional Eulerian frame is not being used for 
the short period which accounts for the fact that pitch attitude dependence 
of the aerodynamics appears in the equations and there is no inertial 
term depending on pitch rate in the normal force equation. 

Z. and M. are the usual downwash lag derivatives. Note that the angle 

Vf 

of attack change of the aircraft is equal to (6 £ - y— ) so that the 

derivatives Z* and M' also include downwash lag effects in addition to 
o p 

the direct effects of pitch rate. 

If we convert these equations of motion into dimensionless form 
consistant with the notation of the rest of this report. 


(1 - Z.) s - Z 
^ w w 


(” 
/ e 


C82') 

- M. s - M 

- «0_ 


f ) 


L w w 





where 

_ '^f 

'^f V£2R 
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z. 

w 

II 

2 

y 

2Y 

aa 

9Cz 

9w. 


m*V 

■f 


Z 




Z 

w 

w 

a 

= n 

C"D 

9w_ 


z. 

0 Vf2R 






M. = M. VR 
w w 


1 ^ ^ 
I* aa 9}<^ 


M 


w 



W 


M. 


Q 


% 



and 


m* = 


M 




— /— I 
2' 2 -^b 


I* = 


" I 
2 b 


C 


Z 


Z 

pTTR^CnR)^ 


M 


M 

p-rrR^Cf^R)^ 


The various derivatives can be expressed as follows: 

1.) The vertical acceleration derivatives (w^) are due 
and arise from the horizontal tail contributions. They 


aa 9w^ 


yV 


de 


iraa ^h\ *^w \ da 


!5l = s c2 a ^ 

aa 9w^ iraa H t w t da 


to downwash lag 
are given by 
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20 The vertical velocity derivatives (w^) arise from the horizontal tail, 
rotors, wing and fuseage. 


aO horizontal tail 


aa 9 w£ 
aa 3wj- 


hT 


S a n - 
Traa H t da^ 


hT 


Traa 




V \ S 


b . ) rotor 


aa 


9C, 


9w, 


aa 


!5i 

9a 

g 


aa 


9C. 


M 


9w, 


2y 


aa 9a 
g 


owing to the equivalence between a gust input a and a vertical velocity 

g 

of the aircraft w- 
t 

c . ) wing 


aa 



w 


^2Sa 

naa w w 


aa 


8C, 


M 


9w , 


w 



traa 


S A c a 
w w w w 


where Z is the distance between the aircraft center of gravity and the 
w 

wing MAC (divided by c^R) 


dO 


fuselage 



aa 9w^ 

2y 


aa 9w^ 
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f ■ 


iraa f 


YV^S 

Traa *f 
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3.) The pitch rate derivatives C0£) arise from the horizontal tail and 


the rotors with a small contribution from the wing 
a.) horizontal tail 





( 

aa 

30f 


iraa H T 

2y 



sr 

a0 

90£ 


TTaa H ^ 

rotors . 

Noting 

the equiv; 


'4 


2y 

aa 

SSf 

T 

aa 3d 

y 

aa 

^0f : 

r 

2y 

aa 3cty 

The wing produces a small 




- "YV 2 S 

aa 

^0f 

w 

TTaa 




- 2 S 

aa 

Hf 

w 

TTaa 




2 2 


de. 


.2 2 . 


The fuselage contribution is neglected. 

4.] The pitch attitude derivatives are directly related to the vertical 
velocity derivatives. 


2y 


aO 30 , 



aa 3w 


Z 

f 


2l % = 2y 

aO 30 _ acr 3wj. 
f f 
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These various derivatives are numerically evaluated using the parameters 
of the aircraft given in Table III to yield the results given in Table IV. 

We now proceed to examine the characteristic equation for the short 
period motion. 

Using the equivalences, = - Zg, = - Mg, and noting that the 
derivative Z^ is negligible compared to 1, the characteristic equation 
for the short period mode is obtained from equation (82') is 




[- 


Z - Ma 
w 0 


" * IVe - “e 




(83) 


Consider the numerical contributions of the derivatives to the stiffness 
and damping terms in this quadratic equation 
c- ZJ . C- Mg.) . (- Z^M.) 

(0.0362) + (0.108) + (- 0.0012) = 0.143 

\ C- ^0^ " C- hV ■ 


(0.00391) + (0.0106) + (0.00039) + (- 0.00121) = 0.0129 


The characteristic roots are 

2 ~ ~ 0.072 ± 0.088 i 

The damping ratio and natural frequency are 



0.63 



0.114 
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Total Aircraft: 

AFT CG (FS 298.2, WL 73.63) ■ 

= 712 (Gross Weight 13000 lbs) 
I* = 81.9 (J^= 12903 slugxft^) 

V - . 844 (300 kts when = 600 

ft/sec) 

Horizontal Tail; 

= .322 
= 4.23 
= 4.13 

- .335 ■ 

da. 

Rotor: 

! 

Wing : 

// = 3 

25„= 1.101 (= 2 e„) 

a = 5-7 

=1*333 

a = .089 


= 105 slugxft^ 

7, = .107 
.w 

y = 3.83 

a,. = 4.58 
w 

I* = = 1.0 

3 3a 


J* = .67 

Fuselage: 

= 1.0355 


= 1-33 

Sj, = .206 

= 3.64 

h = .361 

K = -.344 

ot 

p 

C. = 7.10 

K = -.3 



Table HI. Data Used for the Evaluation of the Derivatives 
Shown in Table IV. 
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Horizontal 

Tail 


Rotors 

(Two) 


Fuselage 


Total 


■ ..J, 

2y 

ao 9 Wy. 

- 1.57 

( * ) 

( * ) 

( * ) 

- 1.57 

Z. 

w 

-.00465 




-.00465 

2y z 

aO 9 

- 1.51 

-.406 X 2 

- 8.64 

- 1.28 

-12.2 

Z 

w' 

-.00446 

-.00240 

-.0255 

-.00379 

-.0362 

2y z 

aa 9 0^ 

6.27 

3.29x2 

-.454 

( * ) 

12.4 


.0185 

.0194 

-.0013 


.0366 

aa 9 dj. 

1.51 

.406 X 2 

8.64 

1.28 

12.2 

✓V 

.00446 

.00240 

.0255 

.00379 

.0362 

- dc 

2>f m 

aa 9w ^ 

2.75 

( * ) 

( * ) 

( * ) 

2.75 

/s. 

M. 

w 

.0335 




.0335 

2y m 

aa 9'.vy 

2.64 

-.175 x2 

-.•382 

- 1.03 

.869 

M 

w 

.0322 

-.00427 

-.00467 

-.0126 

.0106 . 

-X 

Zy' OT 

aa 9 0^ 

- 10.9 

1.05 X 2 

-.0200 

( * ) 

-8.86 

^9 

-.133 

.0256 

-.0002 


-.108 

2, 

aa 9 6^ 

r\ 

M 

e 

- 2.64 
-.0322 

.175 X 2 
.00427 

.382 

.00467 

1.03 

.0126 

-.869 

-.0106 


( * ) : neglected 


Table IV. Derivatives: Results of Numerical Evaluation for the Case 

Given in Table III. 
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This is the quasi-static prediction of the short period characteristic. 

Noting from previous sections that the low frequency rotor flapping mode 
has a frequency approximately equal to 0.18 it can be seen that this 
frequency is reasonably close to the short period motion and would be 
expected therefore that including in some approximate way the coupling of 
this mode to the short period would be the source of the discrepancy in 
damping between the complete model and the quasi-static model described 
previously. 

First consider the relative magnitude of the contributions of the 
various derivatives to the short period motion and in particular the 
rotor contributions. The primary contribution to the damping is Mg and 
next in importance is Z . The influence of ZaM. is small. From Table IV 
it can also be seen that the rotor contribution to these three terms are 
18%, 1.7% and 0.5% respectively. Therefore we can conclude that the important 
contribution of the rotor to the damping term is through Mq, 

In the stiffness term, the most important contribution is Mg. The 

rotor contribution to this tenn is approximately 42% of the total. The 

next term in order of size is Z Ma and from Table IV the rotor contribution 

w 0 

to Z is about 6% and 24% of MA, Therefore it can be assumed that the 

w 0 

primary contribution of the rotor to this term is through Mg. The last 
two terms in the stiffness are relatively small and especially when the 
rotor contributions to them are considered (Table iv) it is ?found that 
they may be neglected. Thus we may conclude that the important rotor 
contributions to the stiffness are through the terms Mg and Mg . 


122 



Now in order to include the rotor dynamics in the short period motion 
we consider the following approach. The rotor contribution to the pitch 
attitude and pitch rate derivatives are expressed as 


aa 


- A Ga (s) 
£ 


where A and- A' are the rotor derivatives when the quasi-static approximation 
is made and Gq and Gq are transfer frunctions which include the effect 
of rotor dynamics. 

Since there is an equivalence such that 


and 


3C ) _ 9( 3 



3 ( ) _ 30 

9a 99 
g 


we can obtain approximations to the dynamic effects by considering the 

rotor hub moment frequency response characteristics to and which 

can be calculated from expressions previously given. 

Figures 48 and .49 show ...the frequency response characteristics of the rotor 

hub moment to sinusoidal inputs in dy and a^. Since many factors 

contribute to these frequency response characteristics .it is difficult 

to obtain simple analytical expressions for the transfer characteristics 

Gg Cs3 and Gq (s) . However, second order approximations can be developed 
£ £ 
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Figure 48. 


Polar Plot: Frequency Response of the Pivot Pitching Moment 
to the Shaft Pitch Rate Input 
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Figure 49. Polar Plot: Frequency Response of the Pivot Pitching Moment 
.to Normal Gust Input. 
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so that 


G- (s) = - C^20I: 

®f + (23 (0.46) (0.20) s + (0.20)^ 


^ ^ , 0.76 + 0.37s + (0.20)^ 

G (s) = ^ 

f s + 2 (0.46) (0.20) s + (0.'20) 


(84) 


The agreement between these approximate transfer functions and the exact 
ones are shown in Figure 50 indicating that a good match is obtained by’ 
second order form up to a frequency of 0.1. 

We can now use root locus techniques to examine the effect of rotor 
dynamics on the short period motion. 

First consider only the rate effect of the rotor written as 




( 0 . 20 )' 


s^ + (2) (0.46) (0.20) s + (0.20)^ 




(M^)r{ 2 

s 


(s) (s + 2 (0.46) (0.20)) 

+ (2) (0.46) (0.20) s + (0.20) 


(85) 


So that the total pitch damping can be written as 


M, - (5) cs 2 C0.46)[0.20)) } 

“ s + (2) (0.46) (0.20) s ♦ (0.20)^ 

is the quasi-static damping. The characteristic equation may be 

written in root locus form as 




(s) (S - z^3 (s + 0.184) 


A [s^ + 2 (0.463(0.20) s + (0.20) 
sp 




= 0 


( 86 ) 
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(a) Pitching Moment to 0^ Input (b) Pitching Moment to 6^ Input 


Figure '50. Polar Plot; Frequency Response of Pitching Moment 
due to Rotor about Aircraft CG to 6^ and 0^ Inputs; 
Comparison with the Approximate- Transfer Functions 


V = .844, Y = 3.83, I* = I*^ = 1.0 
I* = .67, S* = 1.035, = 1,0355, 

= 1.33, h = .361, K « -.344, K 

? > p > p 


C 


a 


-.3 






where is the quasi-static short period motion. 

The root locus shown in Figure 51 shows the effect of the gain indicating 

that the dynamic response of the rotor force to pitch rate and consequently 
coupling between the rotor dynamics and the short period produces a loss 
in damping. 

If we now also include the influence of the rotor dynamics in the 
attitude derivative 1^^ , the characteristic equation can be expressed as 


s j^Cs + 0.184) (s - Z^) + 0.24 (s - 0.78} j 

s^ + 2 (0.46) (0.20) s + 0.20)^ 


This can be placed in root locus form with the ratio (r|— )n constant and the 
results of the influence of these terms on the short period are also shown 


in Figure 51. Using- values from Table IV gives the short period’ d^amics 
shown. The characteristic roots are now 



due to the coupling between the rotor modes and the body modes. The 
approximate results also agree well with the complete dynamic model. 

Unfortunately, the amount of algebra involved precluded finding 
analytical expressions for these transfer functions representing the blade 
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(a) With Rotor Dynamics Reflected in M* 
Only; cf. Characteristic Eqn, CS6)'; 

O : = .0256 



(b) With Rotor Dynamics Reflected in Both Mg 
and MgS cf. Characteristic Eqn. (86); 

G ; = .0256 (@ « .167) 

□ ; Result of Fig. 47 with Complete Rotor 
Dynamics Retained (18 DOF Model) 


Figure 51. Root Locus: Influence of Rotor : Same as Above, but with QS Approximation 

Dynamics on Aircraft Short Period Applied to All Motion DOF except Aircraft 

Mode Characteristics for the Case Longitudinal Rigid Body Motion DOF 

Given in Table III. = 300 kts , 





dynamics. It can be seen that the natural frequency in the denominator 
lies relatively close to the lower flap mode. 

It then appears that the rotor dynamics should be included in 
predicing the short period characteristics of this type of vehicle. It 
is not clear at this time how the flap natural frequency would influence 
this conclusion. Other studies on helicopters have indicated that as the 
fi&p frequency is increased it becomes important to include this coupling 
between blade motion and fuselage motion. 

It is further interesting to note that the trend of the two modes, 
the short period and the lower gimbal mode (B^ ~ I-) ate such that the 
coupling increases with airspeed. It can be seen from Figure 45 that 
the lower gimbal mode frequency decreases with airspeed and the short 
period frequency increases with airspeed and thus there is more significant 
coupling between these two modes as airspeed increases. This trend is 
also supported by the increasing departure of the short period dynamics 
predicted by the quasi-static model compared to the complete model as 
shown in Figure 45, 
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INFLUENCE OF FEEDBACK WITH FUSELAGE FREE 


We now examine the influence of various feedbacks using the complete 
or eighteen degree of freedom model. 

First we consider the effect of body motion feeback such as might be 
employed to improve the stability and control characteristics of the 
aircraft. Figure 52 shows the influence of a combination of pitch rate 
and pitch attitude feedback to the elevator on the dynamics. This feed- 
back of course improves the damping of the phugoid as shown in the figure 
and also increases the short period damping and frequency. The only other 
significant changes occur in the two wing bending modes , where as might be 
expected, as the gain is increased, the frequencies are lowered and the 
spanwise bending mode is destabilized. Infinite gain essentially corresponds 
to the fuselage fixed case and so the trends with increasing gain are essentially 
the reverse of those shown in Figure 46 with respect to freeing the fuselage. 

As indicated in Figure 52 a large value of the feedback gain does produce 
an instability. 

Figures 53 through 59 show the influence of various wing motion feed- 
backs considered earlier for the complete system dynamics. 

In general it may be noted that there is little change from the simpler, 
model with fuselage fixed considered earlier, particularly for any reasonable 
level of gain. Compare fox example Figures-S’3 with Figure 36. It can be 
seen that there is very little change in the influence of wing motion feed- 
backs on the dynamics of the system whether the fuselage is free or fixed. 
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Figure 53. Root Locus: Effect of 9 ^^ Feedback (Proportional); 


18 DOF Model 0 V = .844 (300 kts) 
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Figure 54. Root Locus: Effect of 0^^ Feedback (Proportional); 

18 DOF Model 0 V = .844 (300 kts) 
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Figure 55. Root Locus: Effect of 6^^ Feedback (proportional); 

18 DOF Model @ V = .844 (300 kts) 
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Figure 56. Root Locus: Effect of -»■ 9^^^ Feedback (Proportional); 
18 DOF Model @ V = .844 (300 kts) 
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CONCLUSIONS 


Classical feedback techniques, combined with the introduction of complex 
coordinates are valuable in providing insight into the factors influencing 
the dynamic stability of this complex aeroela.stic system. 

Sin5)lified or lower order dynamic models can be conveniently developed 
using these techniques which provide qualitative insight into the important 
parameters in the problem, however a large number of degrees of freedom 
appear necessary for quantitative accuracy. 

The essential features of the whirl flutter problem, for the range of 
physical parameters studied are shown by considering the proprotor cyclic 
flapping degrees of freedom coupled with wing torsion and spanwise bending. 
The wing chordwise mode is primarily coupled to the collective rotor modes. 
The complexity of the coupling effects is largely associated with the fact 
that the uncoupled wing modes have natural frequencies near or below once 
per revolution. 

Single loop feedbacks of wing motion to cyclic pitch generally appear 
to stabilize one particular wing mode while destabilizing another. 

Adding fuselage degrees of freedom tends to raise the natural fre- 
quencies of the wing modes and to increase the damping of the wing span- 
wise mode while exerting little influence on the damping of the chordwise 
mode. Including fuselage freedom has only a small influence on the effects 
of wing motion feedback. 

Relevant to the stability and control characteristics of the aircraft, 
the flexible modes couple with the short period mode of the vehicle and 
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influence the damping ratio of this mode. Use of the quasi-static assumption 
for the flexible inodes of the vehicle results in an overestimation of the 
short period damping and this coupling tends to increase with flight speed. 
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APPENDIX I 


BASIC THEORY FOR PROPROTOR DYNAMICS IN HIGH INFLOW 
1.1 Introduction 

In this Appendixwe will develop linearized equations of motion for a 
proprotor and its support in axial flow at high inflow ratio, employing a 
relatively simple analytical model. Although the derivation of these equa- 
tions of motion is shown in detail in Ref, 1, it is briefly 
recapitulated in this Appendix to show the basic analytical approach and 
also for use in the main body of the report for a study of the fundamental 
dynamics of the proprotor in high inflow. The derivation of the more 
elaborate model including aircraft longitudinal dynamics is given in Refv 
erences 2 and 13. 

The analytical model employed here consists of a proprotor with three 
blades and the rotor support which is either a pivotted pylon or a canti- 
levered elastic wing. We consider pure out-of-plane ("flap”) and pure in- 
plane ("lag") motion for each blade which is represented by the deflection 
of the blade spar with no torque offset, no droop, and no sweep, retaining 
only lowest modes of their motion. The rotor support motion comprises pylon 
pitch and yaw motion degrees of freedom for the isotropic pylon, and elastic 
bending (in spanwise and chordwise) and torsion motion degrees of freedom 
for the cantilevered wing with only their lowest modes retained. 

For a typical proprotor aircraft operation in airplane cruise mode, the 
rotor induced velocity is negligibly small compared with the forward velocity 
even for the powered flight, owing to the high inflow and low thrust required 
(the ratio of the rotor thrust in airplane cruise to that in hover is inverse- 
ly proportional to the aircraft lift-drag ratio) . If simple momentum 
theory is applied to estimate the ratio of the rotor induced velocity - 
to the forward velocity, it is approximated for high inflow as vfV ~ ^y/(2K^) 
which is typically of order 0.001. When we consider a proprotor in autorota— 
tion where the rotor torque is zero, the rotor produces a negative thrust, 
however, its magnitude is much smaller than that for the powered flight, which 
then makes it more valid to neglect the rotor induced velocity in high inflow. 

From the proprotor dynamics point of view it appears that there is no 
significant difference between the powered flight and the autorotation, and 
since the latter provides a more simplified treatment we assume an autorota- 
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ting proprotor in axial flow at high inflow ratio in the following, and 
consequently the rotor induced velocity is neglected, 

1.2 Rotor Equations of Motion' . 

Rotating Frame 

Consider a model of a proprotor on its support as shown in Fig, 1.1, 
When the rotor support is in motion represented by linear displacements of 
the pivot and angular displacements about it (both assumed small) as defined 
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and the Coriolis inertia coupling of the flap and lag motion is ignored, since 

it is of order yC^/aa and negligible compared with the aerodynamic coupling 

terms of the flap and lag motion which are of order 1 as will be shown later 

(•yW. & yQ \) . In conjunction with this , the Influence of the blade precone is 
t P 

also neglected. The influence of the blade precone becomes important if the 
blade pitch motion"' dynamics , specifically due to the blade pitch control system 
flexibility, is considered; however, its first order effect can be examined by 
introducing the "pitch-lag coupling paratoeter',' which is an approximate treatment 
for a relatively stiff 'control system. 

Nonrotating Frame 

When the number of the blades is equal to or greater than 3, i.e., i7.> 3, 

it is possible to transform the motion variables described in the rotating frame 
to those in the nonrotating frame by applying the Fourier type coordinate trans- 
formation or the "Mulitiblade Coordinate Transformation" (Ref. l) whose advantage 
lies in the simplification it provides in the equations of motion. 

As to the flapping motion, for example, the new degrees of freedom in the 
nonrotating frame are defined as follows: 

B=-^y3 3=^73 cos TZib 

m m 

( 1 . 2 ) 

^ns - I ^ , ^iP/2" I ^ ^ 

m " m 

With these new degrees of freedom, 3^ can be given as 

where the last term On the right hand side appears only if iP is 
summation is from n = 1 to n = (ff - 1) /2 for iP odd and to 
for iP even. 

The relations for the lag motion degrees of freedom are also 
manner. 

In the present case where iP = 3 is assumed, the new degrees of freedom for 
the flap and lag motion in the nonrotating frame are the following: 


(1.3) 

even , and the 
7Z = (JP - 2)/2 
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Flap : 
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Now we will convert the equations (1.1) to those in the nonrotating 
frame. Consider the flap equation. From Eqn, (1.4,1) the time derivatives 
(with respect to the dimensionless time, £2t) of are as follows: 
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Substituting Eqns. (1.4,1) and (1.5) into the flap equation in (1.1) and 
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applying the operations — y (...), rr 7 (..,)cos ip , and— 7 (...)sin ip 
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to the resulting equation with following relations: 
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we obtain the flap equation of motion in the nonrotating frame as follows: 
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where 


I > -^r>T = 77 I COSli) , M = 77 I Sin l|) 

Fo N Z F Flo N z F ^ FIs N ^ F m 

mm m 

( 1 . 8 . 1 ) 

The subscript 6 added to the inertia constants and to the (rotating) natural 

frequency in the coning equation reflects that the different mode shape and 

natural frequency might be used, according to how the blade root is restrained.’ 

That is, if the blade is cantilevered at its root to the hub which is fixed 

to the rotor shaft, then the coning mode and the cyclic flap mode are both the 

elastic bending mode, hence, there is no difference in the mode shapes and 

natural frequencies to be used. However, for the gimballed rotor the cyclic 

flap mode is almost equal to that of an articulated blade which is hinged at 

the center of rotation, i.e., rigid body mode (n. = r) , whereas the coning 

P 

mode is still an elastic bending mode. In such a case the appropriate mode 
shapes and natural frequencies to describe each motion properly might be 
different. 


In a similar manner the lag equations of -motion in the nonrot-ating frame 
may be given as follows: 
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Again the subscript o added to the inertia constants and to the (rotating) 
natural frequency in the collective lag equation indicates that the different 
mode shape and natural frequency might be used. For an autorotating rotor, 
rigid body mode shape (n^ = r) is used for which ~ ^ ^ZO~ 

and now represents the rotor shaft rotational velocity perturbation (with 

reapect to the rotor support) . If there is no shaft speed perturbation 
considered (i.e., perfect rotor rpm governor is assumed) then the mode shape 
and the natural frequency in the collective lag equation of motion are same 
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as those in the cyclic lag equations of motion 


Aerodynamic Flap and Lag Moments 

Now we will find the expressions for the aerodynamic forcing terms on the 
right hand sides of Eqns. (1.7). By making use of the quasi— steady aerodynamics 
and applying the strip theory with the tip loss neglected, the dimensionless 
aerodynamic flap and lag moments on the blade are given as 



"6 ' 



where F and F are the dimensionless section forces resolved into the hub plane, 
and they are related to the section lift and drag (also dimensionless) as 
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(cf. Fig. 1.3). 

In the following evaluation we will use the rigid b^ody mode .shape for 
n„ and n , since the most significant aerodynamic loading occurs near the 
blade tip where the mode shape is very mich like.- that of the rigid body inode. 
Therefore, from the relations given in the above, we have 
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(1.9) 


where the blade chord is assumed constant over the blade span. 

In the trim condition where the rotor is operating in axial flow, the 
velocity components and the angle-of-attack at the section r, are given as 


Uy= r. Up = V + V = y, P = + y2 ^ 

a=8-(j)= 6 - tarT^ (J//v) 


(1.10) 
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for which the moment coefficients in Eqns. (1.9) become 
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( 1 . 11 ) 


where and are the section lift and drag coefficients at trim, 

respectively. 

Now consider small perturbations from the equilibrium. In order to 

evaluate the variations of the aerodynamic flap and lag moments on the blade 

for those perturbations, it is apparently necessary from Eqns. (1.9) to find 

the variations, 6U, 6u, Su^ dC-. and 5(7,, in terms of the motion variables 

T 1, a 

and in the control and gust inputs. 

First, recalling that a rigid body mode shape is assumed for n„ and n , 
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(1.13) 


The variation, 6i7, is obtained from the relation, + Up» 3-s 
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r 6u^+ V 6Up 

■/ q- y2 


(1.14) 


The blade section lift and drag coefficients may be considered as 
functions of the angle-of-attack and Mach number at that section, hence. 




(1.15) 
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where 
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It should he noted that from Eqns. (1.13), (1.14), (1.15) and (1.16) the 
variations, &U^ eventually 'expressed in terms of . the varia- 
tions, . 6u^gj ^^PA^ Therefore, the variational quantities 

in the flap and lag moments on the -blade may be expanded as follows; 
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where the laerodynamic deri'vatives in the above expressions can be found by 
the integration- along the blade span of the integrands which comprise r, V, 

hM > ^do " ^da > ^dl4 

However, the most important and dominant term in each derivative comes from 
the section lift change due to the angle-of-attack change, hence, we retain 
only term in the evaluation of each derivative, using the relation 
/2a = 1/2, although it is a function of the Mach number. 

Thus, the variations of the aerodynamic flap and lag moments may be 
written as follows: 
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and similarly* 
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Comparing Eqns. (1.17) and (1.18), the derivatives are identified as 
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If we use the inflow angle at the effective radius position for the 
evaluation of the integrals involved in the rotor aerodynamic derivatives, 
then the following expressions result: 
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and is usually taken as 3/4. 

Since the trim values of the aerodynamic flap and lag moments and their 
variations from the trim values, as are given in Eqns. (1.11) and (1.17), 
respectively, are derived in the rotating frame, we will convert them to those 
in the nonrotating frame, again applying the Fourier type coordinate trans- 
formation. Noting that the trim values, and and the derivatives, )^®’ 
and )*®» constants and independent of the blade azimuth position - 

(this is owing to the axial flow assumed in the trim condition) , we need the 
transformed expressions only for rela- 

tions given in Eqns. (1.6), (1.5.1) for 6^, and a similar relationship for 
we obtain the following expressions from Eqns- (1.13)! 
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Therefore, the aerodynamic flap and lag moments in the nonrotating frame 
which have appeared on the right hand sides of Eqn, (1.7) are written in 
terms of the motion variables and the control and gust inputs as follows: 
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(1.24) 


Thus the rotor flap and lag equations of motion are now fully described 
in the nonrotating frame by Eqns. (1.7) and (1,24). It should be noted here 
that these linear equations have only constant coefficients as a result of 
the application of the Fourier type coordinate transformation, which is made 
possible by the assumptions on the number of blades (A/ i 3) and the rotor 
operating condition in equilibrium (axial flow) . 


1 .3 Rotor Hub Forces and Moments 

Before developing the equations of motion for the rotor support, we will 
derive the expressions for the rotor hub forces and moments in the nonrotating 
frame which will appear as forcing terms in the support equations. 

Hub Pitch and Yaw Moments 

The rotor hub pitch and yaw moments originate from the bending moment at 

the blade root due to the flapping motion, and with the nonrotating frame 

degrees of freedom they can be written in dimensionless form (normalized by 

^ I, fi^) as follows: 
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(1.25) 
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where a = N^aR^R /iii? = No /'IT. 

Torque, Thrust, H Force and Y Force 

The rotor torque, thrust, H force, and Y force transmitted to the rotor 
hub comprise inertia and aerodynamic contributions. 

First, the inertia contributions to these forces and moment are given 
in coefficient form as follows: 
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where 




The aerodynamic contributions are written as 
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(1.27) 


where F and -F are the dimensionless section aerodynamic forces in the 
X z 

hub plane already introduced before. The expression for the torque above is 
same as given in Eqn. (1.24.4), i.e.. 


aero ~ '“s ‘ 'o’ + V '^^*6 K '^«6 % 

(1.28.1) 
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For. the evaluation of the forces in Eqns. (1.27) we will take a similar 
approach as we did in the derivation of the aerodynamic flap and lag moments. 
That is, expanding the integrals as 
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(1.29) 


and applying the operations of the Fourier- type coordinate transformation, we 
obtain the following: 
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and the derivatives are found as 
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where / ’s and g 's are already defined by Eqns. (1.20). The integrals 

YL YL 

required for the evaluation of the aerodynamic derivatives in Eqns. (1.19) 
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and (1.31) are given as follows: 
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Again, if the effective radius is used for the evaluation of the integrals, 
then the following expressions result; 
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Thus, combining the inertia and aerodynamic contributions, we have 
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1.4 Equations of Motion, for Proprotor and Pivoted Pylon Support 
Support Equations of Motion 

For the pivoted pylon support we consider only its pitch and yaw motion 

about the pivot. The pylon pitch and yaw equations of motion are written in 
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and the forcing terms on the right hand sides above are given in Eqns. (1.25) 

and (1.34) with x , y ^ z and ct degrees of freedom dropped. It should be 
P P P ^ 

noted again that the coefficients of the support equations are all constant. 


Combined Equations of Motion 

Dropping and degrees of freedom also in the rotor equations 

of motion (Eqns. (1.7) & (1.24)), we obtain the combined equations of motion 
for a model considered here in operational form as follows: 
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In the collective flap and lag equations above, and terms are omitted, 

since they only affect the trim values of 3 and t , respectively. It should be 

o o 

noted that the collective flap and lag motion does not couple with the pylon motion 
for the present model, and that the final set of equations of motion has only 
constant coefficients. 


1 . 5 Equations of Motion for Proprotor and Cantilevered Elastic Wing Support 
Support Equations of Motion 

Consider an elastic wing which is cantilevered at its root to a certain 
fixed support (Fig. 1.4). We assume for simplicity that the wing is rectangular 
and has no sweep, no dihedral, and no incidence, and that its motion is repre- 
sented by the elastic bending and torsion of the wing elastic axis. Retaining 
only lowest motion modes, they are expressed as 
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wing vertical displacement 
" chordwise " 

" torsion ” 

" bending mode shape, normalized to at the wing tip 

• " torsion " , " 1 " 

" spanwise bending degree of freedom, positive upward 
" chordwise " , " rearward 

” torsion degree of freedom, positive leading edge up 

If it is further assumed that the pylon is rigidly attached to the wing 
tip such that the pivot is located at the wing tip elastic axis, then the pivot 
linear and angular displacements are expressed in terms of the wing motion as 
follows : 

Xp- , I/p - 0 . Sp 

“x ■ - ’ “a ■ P ' “1/ 

where is the slope of the bending mode shape at the wing tip, and the 

rotor rotational direction is assumed clockwise on the right wing. 
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The wing equations of motion may be written in dimensionless form as 
follows (normalized by ~ : 
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The wing aerodynamic terms on the right hand sides in the above equations 
may be found as 
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where the integrals above are evaluated by using the approximated mode shapes 
for n and £ , l.e., n - u'^ I (v„ i?) and £ - U / 2?) . 


The rotor force and moment contributions are written as follows; 
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where the rotor huh forces and moments are found in Eqns. (1,25) and (1.34) 
with the relations given in Eqns. (1.38) substituted. Again it should be noted 
that the support equations have only constant coefficients. 


Combined Equations of Motion 

Using relations given in Eqns. (1.38), we obtain from Eqns. (1.7), (1.24), 
(1,34), (1.39), (1.40), &. (1.41) the combined equations of motion for the prop- 
rotor and the cantilevered elastic wing support in operational form as follows: 
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where 
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is same as given in Eqn. (1.37.1) 
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Again it should be mentioned that the final, set of the equations of motion for 
the present model has only constant coefficients. 

1.6 Summary 

In this Appendix we have developed linearized equations of motion for a 
simple analytical model of a proprotor and its support in axial flow at high 
inflow ratio. First, rotor blade flap and lag equations of motion which are 
described 4n the rotating frame are Fourier type coordinate transformed to 
those in the nonrotating frame by introducing new motion degrees of freedom 
in that frame. The expressions for the aerodynamic flap and lag moments in- 
voleved in these equations are also obtained in the nonrotating frame. Then, 
rotor hub forces and moments are found in the nonrotating frame, and using 
these hub forces and' moments, rotor support equations of motion are presented 
for two support models, i.e., for a pivotted pylon support and for a canti- 
levered elastic wing support. Finally, the rotor and support equations of 
motion are combined to give a set of linear, second order ordinary differential 
equations with constant coefficients, thus, standard techniques can be applied 
to examine the dynamics of the system, 

It should be mentioned here that we have not considered any coupling bet- . 
ween pitch and flap or lag in the above. This is because their influences can 
be conveniently dealt with by using conventional feedback techniques. 
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ig. 1.1 A Model of Proprotor on Its Support: Definition of Hub 
Forces and Moments, Pylon Linear and Angular Displacements 
and Gust Velocities; Only m-th Blade is shown. 




Fig, \,2 


Blade Flap and Lag Deformations 
(m-.th Blade) 



(Hub Plane) 


• Fig. 1.3 Typical Blade Section; Velocities 
and the Resulting Aerodynamic Forces 
(Looking Inboard) 







Figure 1.5. 
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